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ABSTRACT 


Two  experimental  studies  of  the  natural  convection  characteristics  of 
heated  protrusions  immersed  in  dielectric  liquids  were  conducted.  The 
first  study  used  a  three  by  three  array  of  simulated  20  pin  dual-in-line 
chips  which  were  made  from  aluminum  blocks  with  foil  heaters.  The  second 
set  of  experiments  used  a  three  by  three  array  of  thermal  evaluation 
devices  mounted  on  an  alumina  substrate.  The  devices  were  8.9  mm  square 
chips  which  contained  resistors  and  a  type  of  temperature  sensing 
transistor.  Both  studies  used  an  insulated  Plexiglas  enclosure  with  a  top 
mounted  heat  exchanger  maintained  at  a  constant  10  °C.  Each  array  was 
mounted  on  a  Plexiglas  substrate,  and  spacers  were  used  to  vary  the 
horizontal  distance  from  the  components  to  the  enclosure  wall.  Five 
separate  enclosure  widths  were  used,  with  a  maximum  spacing  of  40  mm. 

The  vertically  oriented  aluminum  blocks  were  tested  with  FC-71  and 
power  levels  ranging  from  0.115  W/chip  to  2.9  W/chip.  The  non-dimensional 
data  obtained  was  used  to  develop  an  empirical  correlation  which  predicts 
Nusselt  number  as  a  function  of  Rayleigh  number  and  enclosure  width.  The 
correlation  was  accurate  to  within  4%  of  the  array  averaged  data,  and  the 
maximum  uncertainty  in  the  Nusselt  number  was  7.4%. 

The  actual  electronic  components  were  tested  with  FC-71,  FC-43,  and 
FC-75.  Power  levels  ranged  from  0.34  W/chip  to  1.48  W/chip.  Again,  the 
data  obtained  was  used  to  develop  a  Nusselt  number  correlation.  In  this 
case  a  better  correlation  of  the  data  was  achieved  using  Grashof  number 
and  enclosure  width.  The  correlation  is  accurate  to  within  2%  of  the 
array  averaged  data.  The  maximum  Nusselt  number  uncertainty  was  4.7%. 
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I .  INTRODUCTION 


A.  THE  ELECTRONICS  COOLING  PROBLEM 

Computer  silicon  chips  continue  to  become  more  powerful  and  smaller 
year  after  year.  However,  the  surface  heat  flux  produced  continues  to 
increase,  and  the  removal  of  this  heat  plays  a  major  role  in  large 
computer  design.  For  example,  Hitachi's  recent  M-880  general  purpose 
computer  uses  a  100  cm'  water  cooled  module  which  dissipates  almost  850  W, 
and  it  is  predicted  that  heat  fluxes  of  nearly  10*’  W/m'  will  be  reached  by 
the  year  2000  (Bar-Cohen,  1991). 

The  challenge  of  removing  this  eleven-fold  increase  in  surface  heat 
flux  is  formidable.  Forced  convection  air  cooling  methods  are  limited  to 
heat  fluxes  of  10^  W/m*  (Bergles,  1991) .  Conduction  cooling  with  a  cold 
plate  and  direct  immersion  dielectric  liquid  cooling  are  the  cooling 
methods  currently  employed  in  large  mainframe  computers.  However,  with 
the  exception  of  the  Cray  series  of  supercomputers,  major  computer 
companies  have  opted  for  various  types  of  generally  complex  cold  plate 
assemblies . 

Extensive  research  into  direct  liquid  cooling  is  being  conducted  for 
a  number  of  reasons.  First,  potential  cooling  schemes  using  dielectric 
fluid  would  be  much  simpler  than  their  conduction/cold  plate  counterparts. 
Second,  the  advantages  and  disadvantages  of  single  phase  natural 
convection,  nucleate  boiling,  and  forced  convection  methods  for  a  wide 
variety  of  available  dielectric  liquids  needs  to  be  studied  to  determine 
their  basic  heat  transfer  characteristics  in  a  simulated  computer  circuit 
board  environment.  Finally,  the  results  of  the  research  should  point  the 
way  for  the  optimum  dielectric  cooling  method,  which  can  then  be  refined 
in  order  to  mal:e  it  a  competitive  means  of  cooling  for  the  next  generation 
of  computers. 


B .  RELATED  RESEARCH 


Park  and  Bergles  (1987)  conducted  natural  convection  experiments  using 
foil  heaters  mounted  both  flush  and  protruding  from  a  circuit  board.  The 
heat  transfer  coefficient  was  measured  for  a  single  flush  heater  with  two 
heater  heights  of  5  mm  and  10  mm.  Heater  widths  varied  from  2  mm  to  70 
mm.  Additional  experiments  used  a  vertical  array  of  two  or  three  heaters 
with  various  distances  between  them.  Combinations  of  heaters  included  two 
or  three  flush  in-line,  two  flush  staggered,  and  two  protruding  in-line. 
Distilled  water  and  R-113  were  the  fluids  used. 

For  the  single  flush  heater,  the  heat  transfer  coefficient  was  found 
to  increase  as  the  heater  width  decreases,  with  this  effect  more 
pronounced  in  the  R-il3.  For  the  in-line  flush  heaters,  the  heat  transfer 
coefficient  was  higher  for  the  bottom  heater,  while  the  opposite  v.’as  true 
for  the  in-line  protruding  heaters. 

Kelleher  et  al.  (1987)  conducted  a  natural  convection  study  of  a  long 
horizontal  protruding  heater  mounted  on  a  vertical  wall  in  a  water  filled 
enclosure.  Heat  exchangers  on  the  bottom  and  top  of  the  enclosure 
maintained  a  constant  temperature.  Results  for  three  separate  heater 
positions  indicated  that  the  Nusselt  number  decreased  as  the  heater 
position  was  raised.  Additionally,  a  flow  visualization  study  revealed 
that  the  flov;  was  divided  into  two  regions.  The  more  active  upper 
buoyancy  driven  region  accounted  for  most  of  the  heat  transfer,  while  the 
more  sluggish  lower  region  was  driven  by  shear  interaction  with  the  upper 
region . 

Joshi  et  al.  (1990)  performed  a  detailed  natural  convection  study  of 
a  vertically  mounted  three  by  three  array  of  heated  protrusions  in  an 
enclosure  filled  with  dielectric  liquid  FC-75.  The  protrusions  were 
horizontally  oriented  rectangular  aluminum  blocks  sized  to  simulate  20  pin 
dual-in-line  (DIP)  packages.  The  top  and  bottom  enclosure  boundaries  were 
heat  exchangers  set  to  maintain  a  constant  temperature  for  various  heater 
power  levels.  Enclosure  width  was  fixed  at  30  mm. 
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Extensive  flow  visualization  revealed  three-dimensional  transport 
which  varied  with  power  level.  As  the  power  level  was  raised,  the  upward 
flow  increased  in  intensity  and  complexity.  Flow  away  from  the  components 
varied  with  time,  and  this  was  confirmed  by  time  history  temperature 
measurement.  Embedded  thermocouples  were  also  used  to  calculate  heat 
transfer  characteristics.  A  correlation  was  developed  to  compute 
component  temperature  from  the  dissipated  power. 

A  follow-on  investigation  by  Joshi  et  al.  (1991)  utilized  a  vertically 
mounted  three  by  three  array  of  vertically  oriented  heated  protrusions  in 
an  enclosure  filled  with  three  different  fluorinert  type  dielectric 
liquids.  Again,  the  upper  and  lower  boundaries  were  constant  temperature 
heat  exchangers.  Enclosure  widths  of  13  mm  and  3  0  mm  were  used  with 
varying  power  levels  in  this  study. 

It  was  found  that  the  top  and  bottom  enclosure  conditions  affected  the 
component  temperatures  to  a  greater  degree  for  the  lower  power  levels. 
The  effect  of  enclosure  width  was  minimal  on  the  resultant  calculations  of 
Nusselt  and  modified  Rayleigh  numbers.  These  non-dimensional  heat 
transfer  characteristics  were  correlated  in  a  similar  manner  to  the 
previous  investigation. 

A  similar  set  of  experiments  using  ethylene  glycol  as  the  fluid  was 
conducted  by  Keyhani  et  al.  (1991).  Five  heated  protrusions  were 
uniformly  spaced  on  a  vertical  wall  inside  an  enclosure  equipped  with  a 
top  mounted  heat  exchanger.  Six  different  enclosure  widths  varying  from 
13.5  mm  to  45  mm  were  tested  at  power  levels  ranging  from  2  W  to  12  W  per 
heater . 

Flow  visualization  of  the  experiments  revealed  primary  flows  along  the 
vertical  walls  separated  by  a  narrow  core  flow  consisting  of  secondary 
flow  cells.  The  heat  transfer  coefficient  of  the  top  and  bottom  heaters 
was  influenced  marl^edly  by  the  power  level  and  enclosure  width.  A  single 
correlation  for  Nusselt  number  versus  modified  Rayleigh  number  was 
developed  that  was  independent  of  heater  location  and  enclosure  width. 


3 


Thesis  experiments  accomplished  by  Aytar  (1991)  and  Matthews  (1991) 
further  established  the  heat  transfer  abilities  of  dielectric  liquids  on 
simulated  electronic  components.  They  both  used  a  three  by  three  array  of 
20  pin  DIP  sized  aluminum  blocks  mounted  vertically  in  an  enclosure  with 
a  top  mounted  heat  exchanger.  Aytar  studied  the  effects  of  enclosure 
width,  power  level,  and  Prandtl  number  on  horizontally  oriented 
protrusions.  Matthews  performed  similar  experiments  with  the  protrusions 
oriented  vertically. 

C .  OBJECTIVES 

The  objectives  of  this  thesis  were  as  follows: 


1.  Complete  collecting  data  on  Matthews'  vertically  oriented 
experimental  setup.  The  component  power  level  and  enclosure  width  were 
varied  using  FC-71  as  the  dielectric  liquid. 

2.  Reduce  the  above  data  into  useful  dimensional  and  non-dimensional 
parameters . 

3.  Utilize  the  non-dimensional  data  in  developing  an  empirical 
correlation  for  the  Nusselt  number  which  takes  into  account  the  effects 
of  enclosure  width  and  Rayleigh  number. 

4.  Using  a  circuit  board  assembly  provided  by  NSWC,  Crane,  rairicate 
an  actual  three  by  three  electronic  component  array  experiment . 

5.  Collect  natural  convection  heat  transfer  data  on  this  circuit  board 
assembly  using  dielectric  liquids  FC-71,  FC-43,  and  FC-75.  The  power 
level  and  enclosure  width  were  varied  for  each  liquid  similar  to  the 
previous  study. 

6.  Reduce  the  above  data  into  useful  dimensional  and  non-dimensional 
parameters . 

7.  Utilize  the  non-dimensional  data  in  developing  a  single  correlation 
for  the  Nusselt  number  which  takes  into  account  the  effects  of  enclosure 
width  and  Grashof  number. 

8.  Based  on  the  above  findings,  make  a  recommendation  for  the  best 
dielectric  liquid  to  use  for  natural  convection  cooling.  Additionally, 
recommend  additional  areas  for  future  research  using  the  selected 
liquid. 
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II.  EXPERIMENTAL  APPARATUS 


A.  TEST  CHAMBER  ASSEMBLY 

The  test  chamber  assembly  consisted  of  the  rectangular  enclosure  and 
heat  exchanger  used  by  Matthews.  The  enclosure  was  constructed  of  25.4  mm 
and  12.7  mm  thick  Plexiglas.  A  3.2  mm  0-ring  sealed  the  boundary  between 
the  enclosure  walls  and  the  bottom  of  the  heat  exchanger.  The  two  were 
assembled  together  with  12  threaded  studs,  washers,  and  nuts.  Plexiglas 
inserts  of  various  widths  were  used  in  order  to  vary  the  spacing  between 
the  circuit  board  assembly  and  the  wall.  Details  of  the  enclosure  are 
shown  in  Figure  1 . 

The  top  mounted  heat  exchanger  was  a  Plexiglas  and  aluminum  single¬ 
pass  type  v;ith  five  rectangular  channels  for  cooling  water  flow.  Heat 
transferred  from  the  dielectric  liquid  was  conducted  to  the  coolant  across 
the  3  mm  thick  aluminum  plate  which  formed  the  bottom  of  the  heat 
exchanger.  Three  thermocouples  embedded  in  the  aluminum  plate  were  used 
for  temperature  measurement.  A  drawing  of  the  heat  exchanger  is  shown  in 
Figure  2. 

B.  SIMULATED  CIRCUIT  BOARD 

The  first  set  of  experiments  utilized  the  same  component  board  used 
by  Matthews.  Nine  aluminum  blocks  arranged  in  a  three-by-three  array  were 
mounted  on  a  12.7  mm  thick  Plexiglas  substrate.  The  numbering  system  for 
the  components  was  unchanged:  bottom  to  top,  right  column  to  left  column. 
Small  10.6  i2  foil  heaters  were  located  between  the  blocks  and  the 
substrate.  Temperatures  were  measured  using  six  thermocouples  per  block. 
Nino  thermocouples  were  also  mounted  on  the  back  of  the  substrate  for 
conduction  loss  calculations.  Details  of  the  heater,  thermocouple,  and 
block  mounting  procedures  are  described  by  Aytar  (1991) .  A  drawing  of  the 
circuit  board  is  shov;n  in  Figure  3. 


5 


25.0  mm 


Section  AA 


Fltrura  2  .  Heat  Exchanger 


blocks: 

7  mm  12.7  mm  — ►  •—  23.9  mm  Ion' 

I  7.9  mm  wide 


C.  NSWC  CIRCUIT  BOARD 

The  remainder  of  the  experiments  used  a  three  by  three  array  of  Texas 
Instruments  thermal  evaluation  devices  which  were  assembled  on  a  50.8  mm 
square  alumina  substrate  board  by  Naval  Surface  Warfare  Center  (NSWC) , 
Crane,  Indiana.  Each  device  is  an  8.9  mm  square  chip  which  contains  four 
resistors  and  a  Temperature  Sensing  Element  (TSE) .  Again,  the  same 
component  numbering  system  is  used  for  consistency.  The  resistors,  when 
connected  in  series,  have  a  resistance  of  approximately  165  ohms.  The  TSE 
is  a  solid  state  temperature  measuring  device,  or  a  type  of  transistor. 
Figure  4  is  a  photograph  of  the  circuit  board,  and  a  schematic  diagrarri  of 
the  chip  internals  is  shown  in  Figure  5. 


Flgurs  4.  NSWC  Circuit  Board 

To  simulate  actual  electronic  chips,  power  is  provided  to  the 
resistors  for  heating  purposes.  Using  a  constant  1  mA  current  source,  the 
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voltage  across  the  transistor  base  to  emitter,  is  measured.  Component 

temperature  is  then  obtained  from  a  previously  plotted  calibration  curve 
of  temperature  versus  ICj..  The  complete  circuit  board  is  mounted  in  the 
center  of  a  Plexiglas  substrate  board  with  essentially  the  same  dimensions 
as  in  Figure  3.  18  thermocouples  were  added  as  follows: 

•  Five  on  the  chip  lid  surfaces.  Located  on  chip  #2,  4,  5,  6  and  9. 

•  Four  on  the  substrate  face.  Located  diagonally  between  chip  #1  and 
5,  #3  and  5,  #"7  and  5,  and  #9  and  5. 

•  Nine  on  the  circuit  board  assembly  back.  Located  directly  behind  the 
chips . 

The  complete  circuit  board  is  shown  in  Figure  6. 


Flgur«  5.  Thermial  Evaluation  Device  Circuitry 


D.  SYSTEM  HARDWARE 

1.  Slmulatsd  Circuit  Board 

Copper -constantan  thermocouples,  0.010  inch  diameter,  were  used 
for  temperature  measurement  on  the  simulated  circuit  board  assembly.  Each 
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foil  heater  was  connected  in  series  to  a  2 . 0  i2  ±  2.5%  resistor.  These 


resistors,  in  turn,  were  connected  in  parallel  to  a  0-100  V,  0-5  A  direct 
current  power  supply.  This  arrangement  allowed  for  a  simple  calculation 
of  heater  power  from  measured  voltages.  Details  of  this  sample 
calculation  can  be  found  in  Matthews'  thesis.  The  thermocouples  and 
heaters  were  connected  to  a  Hewlett-Packard  HP-'497A  Data  Acquisition 
System  (DAS)  .  The  input  and  output  to  the  DAS  was  via  an  HP-9826 
microcomputer.  The  data  channels  were  unchanged  from  Matthews' 
experiments,  and  they  are  repeated  below: 


•  Channels  0-53 

•  Channels  54-56 

•  Channels  57-60 

•  Channel  61 

•  Channels  62-70 

•  Channels  71-75 

•  Channel  76 


Aluminum  block  temperatures 

Heat  exchanger  temperatures 

Back  of  board  assembly  temperatures 

DC  power  supply  voltage 

Foil  heater  voltages 

Back  of  board  assembly  temperatures 

Ambient  temperature 


2.  NSWC  Circuit  Board 

As  before,  identically  sized  copper-constantan  thermocouples  were 
used  v;ith  two  HP-3497A  units.  One  HP-3497A  with  two,  twenty  channel  cards 
was  used  to  measure  thermocouple  and  voltage  data.  The  other  HP-3497A  was 
used  only  as  a  current  source.  It  had  one,  twenty  channel  card  n.odified 
to  supply  the  constant  1  mA  current  to  the  TSEs .  Power  supplied  to  each 
chip  was  easily  calculated  by  multiplying  the  chip  resistor  voltage  by  its 
current.  The  current  was  equal  to  the  precision  resistor  voltage  divided 
by  its  resistance. 

The  chip  resistors  were  powered  from  the  same  direct  current 
power  supply.  However,  all  nine  chips  could  not  be  powered  individually. 
This  was  due  to  the  circuit  board  layout  employed  by  NSWC.  The  center 
column,  chip  #4,  5,  and  6,  could  be  powered  individually,  but  the  chip 
resistors  in  the  two  outer  columns  were  each  wired  in  parallel. 
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Additionally,  due  to  circuit  board  constraints,  TSE  junction  temperatures 
could  not  be  measured  for  the  four  corner  chips.  The  channels  were 
numbered  as  follows: 


Channels 

0-4 

Chip  resistor  voltages 

Channels 

5-9 

Precision  resistor  voltages 

Channel 

10 

DC  power  supply  voltage 

Channels 

11-13 

Heat  exchanger  temperatures 

Channels 

41-44 

Substrate  surface  temperatures 

Channels 

45-49 

Chip  lid  temperatures 

Channel 

50 

Ambient  temperature 

Channels 

51-59 

Circuit  board  assembly  back  t  emper  i.tures 

Channels 

60-64 

Chip  Vp,  voltages 
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III.  EXPERIMEKTAX.  PROCEDURE 


A.  HARDWARE  PREPARATION 

Similar  preparations  were  made  for  both  sets  of  experiments.  Once  the 
component  assembly  was  in  place,  the  proper  sized  spacer  was  inserted  for 
the  particular  run.  Careful  measurements  of  assembly  to  enclosure  wall 
spacing  were  performed  to  ensure  accuracy.  Frequently  a  small  amount  of 
silicone  RTV  was  applied  to  the  corners  of  both  the  component  assembly  and 
the  spacer  to  correct  for  any  small  warpage.  The  remaining  steps  were  as 
f  ol  lov7s  : 


1.  The  enclosure  was  filled  nearly  to  the  top  with  the  proper- 
dielectric  liquid. 

2.  The  heat  exchanger  was  bolted  to  the  enclosure,  with  the  0-ring 
providing  a  seal . 

3.  The  circulating  bath  supply  and  return  lines  were  attached  to  the 
heat  exchanger.  The  bath  unit  was  set  to  6-9  °C  (indicated)  and 
energized.  It  was  found  that  this  temperature  setting  was  required  in 
order  for  the  heat  exchanger  temperature  to  be  10  °C  during  the  runs . 
The  system  was  checked  for  leaks. 

4.  The  foil  heaters  or  chip  resistors  were  energized  with  the  DC  power 
supply . 

5.  The  HP-3497.2i  was  energized.  All  channels  v/ere  scanned  to  insure 
continuity . 

6.  Foam  insulation  v;as  attached  to  the  enclosure  walls  and  bottom. 

7.  Additional  dielectric  liquid  v;as  siphoned  into  the  enclosure  via  the 
vent  hole.  The  air  bubbles  trapped  directly  under  the  heat  exchanger 
were  manipulated  to  the  lead  access  hole  via  a  small  slot  in  the  tops 
of  the  component  assemblies. 


B .  EXPERIMENTAL  PROCEDURE 

After  the  above  preparations  were  completed,  the  proper  voltage  was 
set  on  the  DC  power  supply.  Data  was  taken  after  steady  state  conditions 
were  reached.  For  the  FC-71,  this  took  at  least  eight  hours  from  ambient 
conditions.  .Subsequent  runs  took  two  to  six  hours  to  reach  steady  state, 
depending  on  power  level  and  spacing.  Steady  state  was  achieved  when  two 
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to  three  data  runs  taken  approximately  10  minutes  apart  indicated  a  ±0.1 
°C  random  temperature  difference  between  the  thermocouple  readings.  This 
criterion  was  changed  to  conform  with  Matthews'  work  for  other  dielectric 
liquids.  It  was  ±1  °C  and  ±0.4  °C,  respectively,  for  FC-75  and  FC-43 .  As 
a  rule,  steady  state  conditions  were  achieved  sooner  with  FC-75  and  FC-43. 

For  the  simulated  circuit  board,  data  acquisition  was  accomplished 
with  the  software  program  ACQUIRE.  Calculations  were  then  performed  with 
the  program  CALCDIEL.  Both  of  these  programs  were  originally  written  and 
modified  by  Pamuk,  Benedict,  Torres,  Powell,  Aytar,  and  Matthews.  They 
are  included  in  Appendices  A  and  B. 

For  the  NSWC  circuit  board,  two  new  programs  were  written.  The 
program  ACQ2  was  used  for  data  acquisition,  and  calculations  were 
performed  by  the  program  CALC2 .  Major  improvements  in  the  new  software 
are  as  follows: 

•  Since  the  TSEs  had  to  be  calibrated,  the  thermocouples  were 

calibrated  at  the  same  time.  Calibration  was  performed  using  a 

constant  temperature  oven  and  a  platinum  resistance  thermometer.  The 
equations  for  the  resulting  calibration  curves  were  used  in  the 
programs  to  convert  voltages  to  temperature.  Temperature  uncertainty 
was  calculated  to  be  0.275  °C  for  the  TSEs,  and  it  was  0.3  °C  for  the 
thermocouples . 

•  To  promote  understanding  and  facilitate  future  modifications, 
extensive  documentation  and  explanations  are  included  in  the 
programs,  where  appropriate. 

The  prograrris  ACQ2  and  CALC2  are  included  in  Appendices  C  and  D. 
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IV.  DATA  ANALYSIS 


A.  SIMULATED  CIRCUIT  BOARD 

The  experiments  using  the  FC-71  were  virtually  identical  in  nature  to 
the  experiments  previously  completed  with  FC-75  and  FC-43 .  A  detailed 
explanation  of  the  methodology  used  for  data  analysis  can  be  found  in 
Matthews'  thesis. 

B.  NSWC  CIRCUIT  BOARD 

The  program  CALC2  was  used  to  obtain  the  Nusselt,  Rayleigh,  and 
Grashof  numbers  for  the  NSWC  circuit  board.  These  non-dimensional 
parameters  were  calculated  on  a  single  chip,  horizontal  row,  and  array 
basis  for  three  dielectric  liquids  (FC-75,  FC-43,  and  FC-71)  and  five 
enclosure  spacings.  Various  power  levels  ranging  from  approximately  0.34 
W/chip  to  1.43  W/chip  were  tested.  Assumptions  used  in  the  CALC2  program 
were  as  follows: 


•  The  chip  was  modeled  as  a  square  wafer. 

•  Chip  and  lid  temperatures  of  the  four  corner  chips  were  assumed 
identical  to  the  horizontally  adjacent  chip  in  their  respective  row. 

•  Chip  side  temperature  v.'as  the  average  of  the  TSE  and  the  lid 
temperatures . 

•  Conduction  was  assumed  one  dimensional  from  the  chip  to  the  back  of 
the  Plexiglas  substrate.  The  heat  was  conducted  from  the  chip, 
through  the  alumina  circuit  board  and  a  very  thin  layer  of  silicone 
rubber,  to  the  Plexiglas. 

•  Thermophysical  properties  of  the  above  materials  were  assumed 
constant  at  a  reference  temperature. 

•  Thermophysical  properties  of  the  dielectric  liquids  were  assumed 
constant.  They  were  evaluated  at  T(^i„. 

•  The  temperature  difference  used  for  the  calculation  of  the  heat 
transfer  coefficient  was  area  weighted.  The  lid  accounted  for  about 
55%  of  the  convection  area,  so  the  surface  temperature  was  55%  of  T.^,, 
plus  45%  of  T.j,^,.. 

•  All  contact  resistances  were  assumed  negligible. 
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Calibration  curve  equations  for  the  thermocouples  and  the  TSEs  were 
entered  into  the  program.  All  thermocouple  temperatures  were  obtained 
from  the  following  equation: 

T  (°C)  =  G. 24977483  +  24.896088V  -  0.079219169V' 
where  V  is  the  thermocouple  voltage  in  millivolts.  For  accuracy,  each  TSE 
had  its  own  calibration  curve.  A  representative  equation  for  chip  #2  is 
as  follows: 

T  {°C)  =  577.58074  -  575.54353V\. 

Similar  equations  were  obtained  for  the  other  TSEs.  The  calibration  curve 
for  the  thermocouples  is  shown  in  Figure  7,  and  a  TSE  calibration  curve  is 
shown  in  Figure  8.  The  calibration  data  for  the  five  TSEs  is  included  in 


Figure  7.  Thermocouple  Calibration  Curve 

Appendix  F. 

Due  to  the  wiring  scheme  of  the  NSWC  circuit  board,  the  power  supplied 
to  the  TSEs  had  to  be  calculated  two  different  ways.  The  power  for  the 
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Power  = 


htr  ^ zp 

3Rp 

with  the  correct  (chip  #2  or  8)  substituted  in  the  equation. 

The  heat  loss  by  conduction  involved  several  materials.  The  thermal 
resistance  for  conduction  loss,  R^,  was  represented  as: 


where 


A  =  cross-sectional  area  for  conduction 
L.  =  material  thickness 
k.  =  material  thermal  conductivity 
The  equation  for  calculated  for  each  chip,  was  therefore: 


Q 


loss 


Rc 


where  AT^  was  the  difference  in  temperature  between  the  TSE  derived 
temperature,  T..  .,  and  the  circuit  board  assembly  back  temperature,  Tj,. 

The  net  heat  transferred  from  the  chip  to  the  dielectric  liquid  could 
then  be  calculated  from  the  follow'ing  equation: 

Q.,.  =  Power  - 

The  average  heat  transfer  coefficient,  h,  was  calculated  from: 


where 


h  = 


^net 

A,ot^T 


A,  .  =  total  surface  area  for  convection 

AT  =  area  based  temperature  difference  between  the  chip  surface  and  the 
heat  exchanger,  or  sink 

In  equation  form. 
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AT 


(0.55T,,.  +  0.45T,,.  > 


Density,  p  (kg/m'*) 

FC-7  5:  p  =  (1.825  -  0 . 00246  x  1000 

FC-43:  p  =  (1.913  -  0 . 00218xTfij^^)  x  1000 

FC-71:  p  =  (2.002  -  0 . 00224xrfij„)  x  1000 

Specific  heat,  (J/kg  °C) 

FC-75, 43, 71;  Cp  =  (0.241111  +  3.7037x10'*  x  x  4187 

Kinematic  viscosity,  v  (m*/s) 

FC-75:  V  =  (1.4074  -  2.964x10-2  x  +  3.8018x10-*  x 

-  2.7308x10'®  X  Tjii„  +  8.1679x10'®  x  Tfn„)  x  10  ® 
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FC-A3:  V  =  (8.875  -  0  .  ^1 001  +  1.387xl0'2  x  r|ij„ 

-  2.1469x10'^  X  r|i2„  +  1.3139x10-®  x  TjnJ  x  10'® 


FC-71:  V  =  10-®  X  exp  (6. 8976  -  0.1388x7^2^^  +  1.331x10-^  x 
-  7.041x10“®  +  1.523x10'®  x  Tf^J 


Coefficient  of  thermal  expansion,  p  (1/°C) 


FC-1 5 :  p 


FC-43 :  p 


FC-1 1 :  p 


0.00246 

1.825 

-  0. 00246 xr„2;„ 

0.00218 

1.913 

-  0.00218x7222;:, 

0.00224 

2.002 

-  0.00224x7222;^ 

Nov;  the  various  dimensionless  numbers  which  characterize  the  heat 
transfer  can  be  calculated.  First,  the  ratio  of  thermal  energy  conduction 
to  storage,  or  thermal  diffusivity,  was  found  from; 


Then  the  Prandtl  numiber  could  be  calculated  from: 


Pr  = 


_v 

a 


A  primary  measure  of  convective  heat  transfer,  the  Nusselt  number,  was 
defined  as: 
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where  L  was  the  vertical  length  of  an  individual  chip.  Natural  convection 
effectiveness  is  measured  by  the  Grashof  number,  which  can  be  c,.iculated 
from : 


Gr 


grpL3(r,^  -  W) 
\|2 


Finally,  the  Rayleigh  number  is  defined  as: 

Ra  =  GrPr 
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V.  SIMULATED  CIRCUIT  BOARD  RESULTS 


A .  GENERAL 

The  natural  convection  heat  transfer  characteristics  of  a  vertically 
oriented  array  of  simulated  electronic  components  were  studied  using  the 
dielectric  liquid  FC-71  as  the  coolant.  The  experiments  were  performed  on 
the  same  equipment  used  by  Matthews.  The  enclosure  widths  used,  after 
careful  measurement,  were  determined  to  be  7 ,  9,  16,  28,  and  40  mm.  The 
2  mm  difference  between  these  widths  and  those  reported  by  Matthews  is  due 
to  the  actual  enclosure  width,  with  no  spacers,  being  40  mm  wide  instead 
of  42  mm.  The  same  Plexiglas  spacers  were  used  for  the  FC-71  runs. 
Additionally,  the  same  approximate  power  levels  of  0.115,  0.34,  0.8,  1.3, 
1.7,  2.25,  and  2.9  W/component  were  used  for  the  FC-71  study. 

The  non-dimensional  data  obtained  with  the  FC-71  was  combined  with 
Matthews'  results  with  the  FC-75  and  FC-43 .  An  empirical  correlation  for 
the  Nusselt  number,  Nu,  was  then  derived  in  a  similar  manner  for  the  third 
liquid,  FC-71.  As  defined  in  Matthews'  thesis,  this  correlation  accounted 
for  variations  in  Rayleigh  number  and  chamber  width.  Conspicuously  absent 
is  the  variation  due  to  Prandtl  number,  Pr.  After  reviewing  Matthews' 
results,  it  was  determined  that  the  effect  of  Pr  was  accounted  for  in  the 
Rayleigh  number.  The  correlation  is  of  the  form; 

Nu  =  a  Ra^’ 

where  Nu  is  based  on  the  component  dimension  in  the  direction  of  gravity, 
X  is  a  non-dimensional  enclosure  width,  and  a,  bl ,  and  b2  are  constants. 

B.  DIMENSIONAL  RESULTS 

The  array  average  temperature,  T^.^^  -  is  plotted  against  net 

power,  Qnef  /  iti  Figure  9.  Figure  10  through  Figure  14  show  the  same  plot 
for  all  three  liquids,  with  each  figure  representing  a  different  spacing. 
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The  general  shape  of  the  curves  in  Figure  9  are  identical  to  similar 
graphs  of  FC-43  and  FC-75  data  taken  by  Matthews  (vertical  orientation) 
and  Aytar  (horizontal  orientation) .  However,  the  aluminum  block 
temperatures  are  much  hotter  when  using  FC-71.  Matthews  took  the  maximum 
increase  in  the  array  average  temperatures,  which  occurred  at  2.86  W,  and 
calculated  the  average  of  the  five  spacings.  He  reported  this  temiperature 
to  be  34.9  °C  for  FC-75  and  48.3  “C  for  FC-43.  For  FC-71,  this 
temperature  is  67.2  °C. 

The  maximum  component  temperature  for  FC-71  cooling  was  76.5  °C,  and 
it  occurred  at  a  power  level  of  2.9  W  and  an  enclosure  width  of  7  mm. 
Corresponding  values  for  FC-75  and  FC-43  were  52  °C  and  68  °C, 
respect ively . 

To  be  consistent  with  previous  work,  all  FC-71  component  temperature 
data  was  averaged.  The  temperatures  of  the  three  blocks  on  any  row  were 
averaged  to  facilitate  a  row-by-row  comparison.  For  all  power  levels  and 
spacings,  the  order  of  average  row  temperatures  were  always  top  >  middle 
>  bottom.  This  pattern  indicated  that  the  buoyancy  forces  overcame  the 
viscous  forces  in  the  fluid.  The  boundary  layers  emerging  from  each 
component  were  definitely  affected  by  the  natural  convection  flow  below 
them.  This  pattern  for  FC-71  compares  to  middle  >  top  >  bottom  for  FC-75 
and  top  >  middle  >  bottom  for  FC-43. 

Specific  component  and  row  temperature  extremes  for  the  FC-71  data 
were  as  follows: 

•  Maximum  block  temperature  occurred  on  chip  #3  60%  of  the  time  and 
chip  #6  or  9  34%  of  the  time. 

•  Minimum  block  temperature  occurred  on  chip  #7  97%  of  the  time. 

•  Maximum  temperature  difference  between  the  top  and  middle  rows  was 
1.7  °C.  For  the  middle  and  bottom  rows,  it  was  6.0  °C.  These 
differences  were  noted  at  a  power  level  of  2.9  W  and  a  spacing  of  7 
mm . 
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Flsruro  9.  Array  Average  Temperature  vs.  Net  Power  for  FC-71,  Vertical 
Orientation 
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Figure  10.  Array  Average  Temperature  vs.  Net  Power 
7  rrim  Spacing,  Vertical  Orientation 
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Flgura  12.  Array  Average  Temperature  vs.  Net  Power  for  FC-71,  43,  and  75, 
16  mm  Spacing,  Vertical  Orientation 
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Figure  13.  Array  Average  Temperature  vs.  Net  Power  for  FC-71,  43,  and  75, 
28  rriiri  Spacing,  Vertical  Orientation 
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Fl9ur«  14.  Array  Average  Temperature  vs.  Net  Power  for  FC-^l,  43,  and  75, 
40  mm  Spacing,  Vertical  Orientation 


For  the  0.115  W  power  level,  only  0.5  °C  or  less  separated  all  chip 
temperatures  for  the  7,  9,  and  16  mm  spacings.  This  is  within  the 

uncertainty  of  the  thermocouple  measurement. 

C.  NON-DIMENSIONAL  RESULTS 
1.  General 

The  following  maximums  and  minimums  for  all  power  levels  and 
enclosure  widths  were  noted: 


•  The  maximum  value  of  the  flux  based  Rayleigh  number,  Raf,  was  616.4 

X  10^.  It  occurred  on  chip  #3  at  a  2.9  W  power  level  and  7  mm 

spacing.  The  minimum  Ra,  was  1.417  x  10^  on  chip  #7  at  0.115  W  and  a 
40  mm  spacing. 

•  The  maximum  temperature  based  Rayleigh  number,  Ra ,  was  31.11  x  10’. 
It  occurred  on  chip  #3  under  the  same  conditions  as  the  Ra,  maximum 
listed  above.  The  minimum  Ra  was  190,000  on  chip  #7  at  0.115  W  and 
a  40  mm  spacing. 

•  The  maximum  Nu  was  28.88  on  chip  #4  at  2.9  W  and  a  28  mm  spacing. 

The  minimum  Nu  of  5.40  was  noted  on  chip  #7  at  0.115  W  and  a  7  miri 

spacing . 

•  The  maximum  uncertainty  in  Ra  and  Nu  were  6.90%  and  7.39%, 
respectively.  Both  values  were  calculated  on  chip  #7  at  a  power  of 
0.115  W  and  a  spacing  of  9  mm. 

Both  Aytar  and  Matthews  plotted  the  array  averaged  Nu  as  a 
function  of  either  the  array  averaged  Ra,  or  Ra.  When  the  data  taken  at 
0.115  W  was  omitted,  it  was  found  that  these  plots  were  straight  lines, 
independent  of  enclosure  width.  This  was  done  since  the  uncertainty  was 
highest  at  this  power  level,  and  the  resulting  chip  temperatures  were  less 
than  or  equal  to  the  ambient  temperature.  Figure  15  is  a  similar  plot  for 
FC-71.  For  comparison,  Nu  versus  Ra  for  all  three  dielectric  liquids  have 
been  plotted  for  each  enclosure  width  in  Figure  16  through  Figure  20.  In 
these  figures,  it  is  noted  that  the  slopes  of  any  particular  dielectric 

liquid  are  virtually  identical  in  nature. 

2.  Effect  of  Rayleigh  Number 

Matthews  and  Aytar  both  observed  a  linear  relationship  when  log 
Nu  was  plotted  against  log  Ra .  This  relationship  was  of  the  following 
form : 
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Rayeigh  Number 


Flgura  15.  Nu  vs.  Ra  for  FC-71,  Array  Averaged,  Vertical  Orientation,  and 
all  Enclosure  Widths 
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Flgur«  16.  Nu  vs.  Ra  for  FC-71,  43,  and  75,  Array  Averaged,  Vertical 

Orientation,  and  7  mm  Enclosure  Width 
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Figure  17.  IJu  vs.  Pa  fcr  FC-ti,  43,  and  75,  Array  .Averaged,  Vertical 
Orientation,  and  ?  rriri  Enclosure  Width 
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Nu  vs.  Ra  for  FC-71,  43,  and  75,  Array  Averaged,  Vertica 

and  16  min  Enclosure  Width 
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Nu  =  C;Ra“' 

where  c-  and  bl  are  constants.  Both  row  averaged  and  array  averaged  values 
were  calculated  and  plotted  for  each  spacing.  The  curve  fit  software 
TABLECURVE  (1990)  was  used  to  find  the  coefficients  of  the  above  equation, 
and  the  software  package  SIGMAPLOT  (1989)  was  used  to  produce  the  graphs. 
Again,  for  consistency  and  correlation  accuracy,  the  data  taken  at  0.115 
W  was  omitted. 

Figure  21  through  Figure  25  are  the  FC-71  data  array  averaged  Nu 
versus  Ra  plots  for  the  five  spacings.  The  corresponding  curve  fit 
equation  is  also  included.  For  completeness.  Figure  26  through  40  are 
similar  plots  for  the  spacings  using  row  averaged  data.  The  value  of  the 
constant  bl  was  found  to  vary  from  0.225  to  0.280  for  the  array  averaged 
data.  The  average  value  of  bl  was  0.249.  Corresponding  values  of  0.381 
and  0.371  were  reported  by  Matthews  for  FC-75  and  FC-43,  respectively. 
The  constant  C;  varied  from  0.255  to  0.514  for  the  same  array  averaged 

data,  and  the  average  value  of  c,  was  0.383. 

3.  Effect  of  Enclosure  Width 

Enclosure  width  effects  on  Nu  were  accounted  for  in  the  following 

equat ion : 

Nu  =  c-.X'^' 


where  b2  is  a  constant  and  c^  is  the  Ra  dependence.  X  is  the  non- 
dimensional  enclosure  width,  which  is  the  ratio  of  the  actual  spacing  to 
the  maximum  spacing  of  4  0  mm.  The  constant  C;.  can  be  represented  by  : 

c^  =  C;Ra^‘ 


The  constant  b2  was  derived  from  a  plot  of  Nu  vs .  X  for  the  five  spacings, 
which  is  shown  in  Figure  41.  The  values  for  Nu  were  taken  from  the  curve 
fit  equations  for  each  spacing's  array  averaged  data  using  an  average  Ra 
of  4  X  10''.  The  resulting  value  for  the  exponent  b2  was  0.165. 

Combining  the  above  results  gives  the  following  general 
correlation  for  FC-71: 

Nu  =  0.383Ra'' ‘‘'''X' 
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RA/Ll-'_>H  NUMBEk 


Flgurs  21.  Nu  vs.  Ra  for  FC-71,  Array  Averaged,  Vertical  Orientation 
and  7  mill  Spacing 


Pigurs  23.  Nu  vs.  Ra  for  FC-71,  Array  Averaged,  Vertical  Orientation 
16  rriTi  Spacing 
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Figure  24.  Mu  vs.  Pa  for  FC-71,  Array  .Averaged,  Vertical  Orientation, 
28  mm  Spacing 


Figure  26.  Nu  vs.  P.a  for  FC-71,  Bottom  Rov.-  Averaged,  Vertical 
Orientation,  and  7  miri  Spacing 
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Figure  27.  Nu  vs.  Ra  for  FC-71,  Middle  Row  Averaged,  Vertical 
Orientation,  and  7  mrr,  Spacing 
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Flsrure  28.  Nu  vs.  Ra  for  FC-71,  Top  Row  Averaged,  Vertical  Orientation 
and  7  mm  Spacing 
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Figure  29.  Nu  vs.  Ra  for  FC-71,  Bottom.  Row  Averaged,  Vertical 
Orientation,  and  9  mm  Spacing 
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Figure  30.  Nu  vs.  Ra  for  FC-71,  Middle  Rov;  Averaged,  Vertical 
Orientation,  and  9  mrv.  Spacing 


Figure  31.  Nu  vs.  Ra  for  FC-71,  Top  Row  Averaged,  Vertic 
and  9  miTi  Spacing 
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Flgura  32.  Nu  vs.  Ra  for  FC-71,  Bottom  Rov;  Averaged,  Vertical 
Orientat.Lon,  and  16  mrr.  Spacing 
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Figur«  33.  Nu  vs.  Ra  for  FC-71,  Middle  Row  Averaged,  Vertical 
Orientation,  and  16  miri  Spacing 


Plgur«  34.  Nu  vs.  Ra  for  FC-71,  Top  Row  Averaged,  Vertical  Orientation, 
and  16  nuri  Spacing 


Figure  35.  Nu  vs.  Ra  for  FC-71,  Bottom  Rov.-  Averaged,  Vertical 
Orientation,  and  28  mm  Spacing 
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Fl^r«  36.  Nu  vs.  Ra  for  FC-71,  Middle  Roa-  Averaged,  Vertical 
Orientation,  and  28  irm'i  Sparing 


Plgur*  37.  Nu  vs.  Ra  for  FC-71,  Top  Row  Averaged,  Vertical  Orientation, 
and  28  mrr,  Spacing 
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Ispira  41.  Nu  vs.  X  for  FC-71, 
11  Enclosure  Widths 


Array  Averaged,  Vertical  Orientation, 


It  is  valid  over  the  ranges; 

9  X  10^  <  Ra  <  2  X  10" 

0.175  <  X  <  1.0 

However,  this  correlation  was  found  to  be  accurate  only  to  within  11%  of 
the  array  averaged  curve  fit  equations.  In  order  to  achieve  improved 
accuracy,  representative  values  of  Ra  were  chosen  for  the  general 
correlation  and  each  spacing's  correlation.  The  Nu  number  results  were 
then  compared,  and  a  trial  and  error  approach  was  used  to  come  up  with  a 
better  value  for  c, .  When  the  value  for  Cj  was  changed  from  0.3  83  to 
0.435,  the  accuracy  improved  to  4%.  Therefore,  a  better  correlation  for 
FC-71  is: 


Nu  =  0.435Ra--'‘'X°'“ 

It  is  valid  over  the  same  ranges  as  listed  above. 
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VI.  NSWC  CIRCUIT  BOARD  RESULTS 


A.  GENERAL 

The  next  logical  step  in  the  dielectric  liquid  natural  convection 
studies  was  to  replace  simulated  electronic  components  with  the  actual 
electronic  devices.  A  three  by  three  array  of  8.9  mm  square  thermal 
evaluation  devices  were  tested  on  a  Plexiglas  circuit  board  assembly  of 
the  same  dimensions  as  the  previous  experiment.  All  three  dieleccric 
liquids  were  tested.  Additionally,  the  same  enclosure,  spacers,  and  data 
acquisition  equipment  were  used  again.  Due  to  manufacturing  tolerances, 
the  enclosure  widths  this  time  were  8,  10,  16,  28,  and  40  mm.  Power 
levels  were  chosen  to  match  those  of  the  previous  studies,  but  they  were 
limited  to  the  125  °C  maximum  temperature  of  the  temperature  sensing 
elements  (TSEs)  inside  each  chip.  Therefore,  TSE  temperatures  were 
purposely  limited  to  about  100  °C  to  meet  this  criterion  in  conjunction 
with  the  TSE's  calibration  raiige.  For  FC-43  and  FC-75,  power  levels  of 
0.34,  0.57,  0.8,  1.3,  and  1.48  W/chip  were  chosen.  For  the  typically 
hotter  FC-71,  the  power  levels  were  0.34,  0.57,  0.8,  and  1  0  W/chip.  Note 
that  the  0.115  W  power  level  was  completely  eliminated  due  to  lessons 
learned  from  the  previous  experiments. 

A  total  of  7  0  data  runs  were  recorded  and  analyzed  for  the  three 
dielectric  liquids.  The  analysis  of  the  dimensional  results  is  similar  to 
Matthews'  and  Aytar's  work,  but  the  non-dimensional  analysis  is  somewhat 
different.  After  the  FC-75  and  FC-43  data  runs  were  complete,  it  was 
realized  that  plots  of  Nu  versus  Grashof  number,  Gr,  were  linear, 
independent  of  Pr  or  power  level.  Therefore,  CALC2  was  altered  to  compute 
Gr  in  parallel  with  Ra .  Average  Pr  was  also  calculated.  The  data,  stored 
on  floppy  disk  with  the  ACQ2  program,  was  then  re-run  with  the  modified 
CALC2  program.  The  extensive  non-dimensional  data  thus  obtained  was  used 
to  produce  a  correlation  i.or  Nu .  The  correlation  is  of  the  following  form: 
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Nu  =  aGr’^'X^^' 

where  Nu  is  based  on  the  chip  dimension  in  the  direction  of  gravity.  X  is 
again  the  non-dimensional  enclosure  width,  and  a,  bl,  and  b2  are 
constants . 

B.  DIMENSIONAL  RESULTS 

1.  PC-71 

As  before,  was  plotted  against  for  all  spacings. 
As  it  can  be  see’-  in  Figure  42,  the  daua  is  linear  for  all  spacings.  The 
maximum  chip  temperature  of  ‘5.4  °C  occurred  in  chip  #6  at  a  power  level 
of  1.0  VJ  and  a  spacing  of  8  mm.  For  comparison  purposes,  the  chip 
temperature  data  was  also  row  averaged  like  the  previous  experiments. 
Hov.’ever,  recall  that  temperature  data  for  the  top  and  bottom  rows  or  TSEs 
WciS  based  on  a  single  chip  each. 

For  all  spacings  and  power  levels,  the  order  of  average  row- 
temperatures  were  always  top  >  middle  >  bottom..  As  in  the  previous  stud'/, 
this  pattern  indicated  the  dominance  of  the  buoyant  forces  over  the 
viscous  forces  in  the  fluid.  The  maximumi  temperature  difrerence  betv;een 
the  top  and  m.iddle  rows  was  1.4  °C.  For  the  middle  and  bottom  row’s,  it 
was  5.4  °C.  These  differences  were  noted  at  a  power  level  of  1.0  W  and  a 
spacing  of  8  mm. 

The  thermocouples  installed  on  the  a2  jmina  s  'bstrate  surface  and 
the  back  of  the  Plexiglas  circuit  board  assembly  also  indicated  a 
temperature  gradient.  The  upper  two  substrate  temperatures  were  within 
0.5  °C  of  each  other;  the  same  could  be  said  for  the  lower  two  temperature 
readings.  The  upper  pair  of  readings  were  always  higher  than  the  two 
substrate  readings  below.  Similarly,  v/hen  the  nine  backside  temperatures 

were  grouped  by  row,  the  magnitude  order  was  always  top  >  middle  >  bottom. 

2.  PC-43 

Results  for  FC-43  were  very  similar  to  FC-71.  Figure  43  shows 
the  T„,,g  -  Tj  versus  Q^,,,  plots  tor  all  spacings.  It  is  noted  that  the 
slope  of  the  FC-43  data  is  smaller  that  its  FC-71  counterpart.  This 
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Piguro  42.  Array  Average  Temperature  vs.  Net  Power  for  F' 
Circuit  Board 
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Flgur«  43.  Array  Average  Temperature  vs.  .Net  Power  for  FC-43 
Circuit  Board 


indicates  that  FC-43  is  a  more  efficient  heat  transfer  medium  for  the 
chips  at  a  given  power  level. 

The  average  row  temperature  order  remained  the  same  as  the 
previous  fluid.  The  maximum  chip  temperature  of  101.9  °C  was  noted  on 
chip  #6  at  a  spacing  of  8  mm  and  a  power  level  of  1.48  W.  Maximum 
temperature  differences  were  2.1  °C  between  the  top  and  middle  rows  and 
4.9  °C  between  the  middle  and  bottom  rows.  As  before,  these  differences 
occurred  at  the  highest  power  level,  1.48  W,  and  minimum  spacing,  8  mm. 
Finally,  the  remaining  substrate  and  circuit  board  assembly  temperatures 

exhibited  the  same  trends  as  described  for  the  FC-71. 

3.  FC-75 

The  third  fluid,  FC-75,  exhibited  almost  identical  results  to  FC- 
43.  In  Figure  44,  the  T,,^^  -  Ts,,,|,  versus  plots  are  again  linear,  and 

the  slope  of  the  FC-75  data  is  the  smallest  of  the  three  liquids  studied. 
The  maximum  chip  temperatures  achieved  was  only  92.3  °C.  As  expected,  it 
occurred  on  chip  #6  under  maximum  power  and  minimum  spacing  conditions. 

As  previously  noted,  the  average  row  temperature  data  was  again 
always  top  >  middle  >  bottom.  However,  the  maximum  differences  between 
rows  was  small.  Only  0.4  °C  separated  the  top  and  middle  row’s,  and  3.3  °C 
was  Whe  difference  between  the  middle  and  bottom  rows.  These  row  averaged 
temperature  differences  were  noted  at  a  spacing  of  8  mm  and  a  power  level 
of  1.48  W.  The  established  pattern  of  top  >  middle  >  bottom  was  again 
observed  fof*  row  wise  substrate  surface  and  circuit  board  assembly 
temperatures . 

4.  FC-71,  FC-43,  and  FC-75  as  a  Group 

For  co>toarison  purposes.  Figure  45  through  Figure  49  are  plots 
of  T„^,j,  -  T^j^j  versus  for  all  three  fluids.  This  means  of  presenting 

the  data  highlights  the  following  conclusions: 

•  For  a  given  power  level,  the  dielectric  fluid  FC-75  convects  heat 
away  from  the  chips  more  efficiently  than  either  FC-43  or  FC-71. 
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Tavg  --  Tsink  (C) 


Qnet  (W) 


Figure  44.  Array  Average  Temperature  vs.  Net  Power  for  FC-75, 
Circuit  Board 


Qnet  (W) 


Figure  45.  Array  Average  Temperature  vs.  Net  Power  for  FC-71,  43,  and  75 
8  mm  Spacing,  NSV.'C  Circuit  Board 


Qnet  (W) 


Figure  46.  Array  Average  Temperature  vs.  Net  Power  for  FC-71,  43,  and  75, 
10  miTi  Spacing,  NSWC  Circuit  Board 
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Figure  47.  Array  Average  Temperature  vs.  Net  Power  for  FC-~1,  43,  and  75, 
16  miTi  Spacing,  HSWC  Circuit  Board 
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Figure  48.  Array  Average  Temperat 
28  miTi  Spacina,  NSWC  Circuit  Board 


Tavg  —  Tsink 


Qnet  (W) 


Figure  49.  Array  Average  Temperature  vs.  Net  Power  for  FC-71,  43,  and  75 
40  mm  Spacing,  NSWC  Circuit  Board 


•  Operation  with  FC-71  leads  to  the  highest  chip  temperatures.  This 
would  be  a  major  disadvantage  in  the  selection  of  the  best  dielectric 
liquid  for  immersion  cooling. 

C .  NON-DIMENSIONAL  RESULTS 

1 .  General 

The  analysis  of  the  non-dimensional  data  for  the  three  dielectric 
liquids  was  grouped  together.  This  was  done  because  each  liquid  behaved 
similarly  for  a  given  power  level  or  spacing.  The  maximum  Nu  and  Gr  were 
always  noted  at  the  highest  power  level  for  a  given  liquid,  which  was  1.0 
W  for  the  FC-71  and  1.48  W  for  the  FC-43  or  FC-75.  Similarly,  the  minimum 
Nu  and  Gr  always  occurred  at  the  minimum  power  level  of  0.34  W  for  all 
liquids.  Additionally,  the  maximum  uncertainty  in  Nu  and  Gr  was  observed 
at  this  minimum  power  level  and  a  minimum  spacing  of  8  mm. 

Specif ically,  the  following  extremes  were  noted  as  described 

above : 

•  Maximum  Nu  at  a  spacing  of  40  mm:  16.37  for  FC-71,  31.04  for  FC-43, 
and  41.44  for  FC-75. 

•  Maximum  Gr  at  a  spacing  of  8  mm:  3450  for  FC-71,  2.478  x  10'  for  FC- 
43,  and  1.224  x  10'  for  FC-75. 

•  Minimum  Nu  at  a  spacing  of  8  mm:  9.81  for  FC-71,  17.19  for  FC-43, 
and  21.32  for  FC-75. 

•  Minimum  Gr  at  a  spacing  of  40  mm:  59  for  FC-71,  1.030  x  10'*  for  FC- 
43,  and  1.112  x  10'  for  FC-75. 

•  Maximum  uncertainties  in  Nu  and  Gr,  respectively:  3.4%  and  2.0%  fcr 
FC-71,  4.2%  and  3.3%  for  FC-43,  and  4.7%  and  3.8%  for  FC-75. 

2.  Effect  of  Grashof  Number 

As  previously  stated,  the  final  goal  of  the  non-dimensional 
analysis  was  to  produce  an  empirical  correlation  for  Nu .  Figure  50 
through  Figure  54  are  plots  of  Nu  versus  Gr  for  each  individual  spacing. 
When  all  three  dielectric  liquids  are  plotted  together  as  shown,  a  linear 
relationship  is  realized. 

This  relationship  was  assumed  to  be  of  the  form: 

Nu  =  CiGr*" 
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d3akNnN  ii3SsnN 


Figure  50.  Nu  vs.  Gr  for  FC-71,  43,  and  75,  Array  Averaged,  NSWC  Circuit 
Board,  and  8  mit  Enclosure  Width 
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Plgura  52.  Nu  vs.  Gr  for  FC-71,  43,  and  75,  Array  Averaged,  NSWC  Circuit 
Board,  and  16  mrri  Enclosure  Width 
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Figure  54.  Nu  vs.  Gr  for  FC-71,  43,  and  75,  Array  Averaged,  NSVrc  Circui 
Board,  and  4C  mrri  Enclosure  V.’idth 


where  c.  and  bl  are  constants.  Again,  TABLECURVE  was  used  to  find  the 
constants  in  the  above  equation.  The  constant  c-  varied  frorri  5.30  to  6.04, 
and  the  average  was  equal  to  5.61.  The  constant  bl  was  found  to  vary  from 

0.121  to  0.141,  with  the  average  being  equal  to  0.133. 

3.  Effect  of  Enclosure  Width 

Enclosure  width  effects  were  accounted  for  by  the  following 

equat ion : 

Nu  =  CjX*^^ 

where  b2  is  a  constant  and  c,  is  the  Gr  dependence.  X  is  again  the  non- 
dimensional  enclosure  width.  Simdlar  to  the  previous  stud-j-,  the  constant 
c_  can  be  represented  by: 

C;  =  c.Gr'^’ 

The  constant  b2  was  derived  from  a  plot  of  Nu  versus  X  for  the  five 
spacings,  v;hich  is  shown  in  Figure  55.  The  values  for  Nu  were  derived 
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Figure  55.  Nu  vs.  X  for  FC-71,  43,  and  75,  Array  Averaged,  NSWC  Circuit 
Board,  and  all  Enclosure  Widths 
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from  the  curve  fit  equations  for  each  spacing's  data  using  an  average  Gr 
of  9  X  10’.  The  resulting  value  for  the  exponent  b2  was  0.154. 

Combining  the  above  results  gives  the  following  general 
correlation  for  the  three  dielectric  liquids: 

Nu  =  5.61Gr°'”X“‘^'‘ 

It  is  valid  over  the  ranges: 

1  X  10^  <  Gr  <  8  X  10‘ 

0.20  <  X  <  1.0 

Five  representative  Grashof  numbers  per  spacing  were  selected  to  determine 
values  for  Nu  from  the  respective  curve  fit  equations.  These  results  were 
then  compared  to  the  Nusselt  number  produced  from  the  correlation.  The 
average  difference  in  Nu  for  the  25  data  points  was  12.3%.  Similar  to  the 
previous  study,  a  trial  and  error  approach  was  used  to  improve  the 
accuracy  of  the  general  correlation.  When  the  value  for  Cj  was  changed 
from  5.61  to  6.40,  the  agreement  improved  to  less  that  2%.  Therefore,  a 
refined  correlation  for  the  three  dielectric  liquids  is: 

Nu  =  6.40Gr'''’X''  =  ''' 

It  is  valid  over  the  ranges: 

1  X  10'‘  <  Gr  <  8  X  10" 

0.20  <  X  <  1.0 
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VII.  CONCLUSIONS 


Two  studies  of  the  natural  convection  heat  transfer  of  heated 
protrusions  immersed  in  dielectric  liquids  were  conducted.  The  first 
study  used  a  three  by  three  array  of  computer  chip  sized  aluminum  blocks 
immersed  in  FC-71.  The  other  study  used  a  three  by  three  array  of  8.9  mm 
square  thermal  evaluation  devices.  Three  fluids,  FC-71,  FC-43,  and  FC-75 
were  evaluated.  Both  studies  used  an  insulated  Plexiglas  enclosure  with 
a  top  mounted  heat  exchanger.  Spacers  were  used  to  vary  the  enclosure 
width,  and  the  maximum  spacing  was  40  mm.  Conclusions  from  the  two 
studies  are  as  follows: 


1.  For  the  first  study,  an  empirical  correlation  for  Nusselt  number  was 
developed.  It  took  into  account  variations  in  Rayleigh  number  and  non- 
dimensional  enclosure  width,  X.  The  correlation  was  based  on  array 
averaged  data.  It  is  listed  below: 

Nu  =  0.435Ra'-‘‘'-X°='^ 

9  X  10’  <  Ra  <  2  X  10" 

0.175  <  X  <  1.0 

The  maximum  uncertainty  in  the  Nusselt  number  was  7.4%,  and  the 
correlation  was  accurate  to  within  4%  of  the  array  averaged  data. 

2.  When  the  FC-71  data  was  combined  with  Matthews'  FC-43  and  FC-75 
data,  several  generalizations  could  be  made.  For  array  averaged  plots 
of  Tj.,,  -  T,..,^  versus  Qn*,. ,  the  order  of  temperatures  for  the  three  fluids 
was  FC-71  >  FC-43  >  FC-75.  Additionally,  the  temperatures  increased  as 
power  level  increased  or  enclosure  width  was  decreased.  For  log- log 
plots  of  Nu  versus  Ra,  each  liquid  exhibited  a  linear  relationship.  The 
FC-71  data  had  the  highest  Nu  and  Ra  numbers  for  each  spacing. 

3.  For  the  electronic  chip  study,  a  general  correlation  was  developed 
for  Nusselt  number  from  the  combined  data  of  the  three  dielectric 
fluids.  This  correlation  took  into  account  variations  in  Grashof  number 
and  non-dimensional  enclosure  width,  X.  As  before,  it  was  based  on 
array  averaged  data.  The  correlation  is  as  follows: 

Nu  =  6 .40Gr''  ‘^^X'’-'-'' 

1  X  10^  <  Gr  <  8  X  10" 

0.20  <  X  <  1.0 


The  maximum  uncertainty  in  the  Nusselt  number  was  4.7%,  and  the 
correlation  was  accurate  to  within  2%  of  the  array  averaged  data. 


4.  Plots  of  versus  for  the  three  liquids  again  showed 
that  the  order  of  temperatures  was  FC-71  >  FC-43  >  FC-75  for  each 
spacing.  For  all  three  liquids,  the  order  of  the  averaged  temperature 
data  was  always  top  >  middle  >  bottom.  This  same  order  was  also 
exhibited  by  the  row  averaged  Plexiglas  substrate  back  and  alumina 
substrate  surface  temperatures. 

5.  Overall,  the  best  liquid  for  natural  convection  heat  transfer  was 
FC-75.  The  best  liquid  is  defined  as  the  one  which  produced  the  lowest 
component  temperatures  for  a  given  power  level  or  spacing.  Lower  chip 
temperatures  equate  to  longer  chip  lives. 
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VIII.  RECOMMENDATIONS 


The  following  recoininendations  are  made  for  further  research: 


1.  Manufacture  additional  NSWC  circuit  boards  that  have  every  chip 
wired  individually.  Use  these  boards  to  assemble  a  large  array 
equivalent  in  size  to  a  typical  mainframe  computer  circuit  board. 

2.  Test  the  above  array  in  a  suitably  sized  enclosure  filled  with  FC- 
75.  Produce  an  empirical  correlation  for  the  Nusselt  number. 

3  .  In  parallel  with  the  experimental  work,  produce  computer  models 
which  can  also  be  used  to  predict  heat  transfer  characteristics. 
Compare  these  results  with  the  experimental  results  and  modify  the 
programs  as  necessary. 


74 


APPENDIX  A.  COMPUTER  PROGRAM  ACQUIRE 
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1670  PRINT  USING  " 1 0X  , " "TOP : " "  .24X  .3( DD. DD  ,£X  1 " ; T{ 54  )  ,T  ; 55 >  ,T ( 56  ) 

1680  PRINT 
1690  PRINT 
1700  I 

1710  PRINT  USING  "5X,""BACK  PLANE  TEMPERATURES  ARE  : . 

1720  PRINT 

1730  PRINT  USING  " 1 0X . " "T( 57  ) : " “ .2X  ,DD. DD" ; T( 57  ) 

1740  PRINT  USING  " 1 0X , " "T ( 58 > : " " ,2X .DD . DD " ; T(  58 ) 

1750  PRINT  USING  "  10X  , '  "^(59  ):  “•■  ,2X  .DD.DD"  :T(SS  ; 

1760  PRINT  USING  " 1 0X , " "T) 60  ) : " "  ,2X  ,DD . DD" ; T( 60  ) 

1770  PRINT  USING  " 1 0X , " "T(  71  ) : ' "  ,2X  ,DD. DD" ;T' 71  ) 

1780  PRINT  USING  " 1 0X  . " " T( 72  ) : " "  ,2X  .DO . DD" : T( 72 > 

1790  PRINT  USING  " 1 0X , " "T( 73  ) : " " ,2X  ,DD . DD" i T( 73  ) 

1800  PRINT  USING  " 1 0X , " " T ( 74 ) : " " ,2X ,DD .DD" ! T( 74 ) 

1810  PRINT  USING  "  1  0X , " "T(  75 ) : “ " .2X .DD.DD" :T( 75  ) 

1820  BEEP 

1630  PRINTER  IS  1 

1840  I 

1850  ASSIGN  eC'le  TO  • 

Ifi60  END 
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40 
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280 

290 
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460 
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490 
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PROGRAM  CalcDiel 


MODIFIED  BY  LOOP  R.  THOMPSON  1/13  ANO  4/15/92. 
FROM  ORIGINALS  OF  PAMUK  ,  BENEDICT,  TORRES, 
AYTAR  AND  MATTHEUS. 


THIS  PROGRAM  ANALYZES  THE  DATA  READ  FROM  • 
A  DATA  FILE  DESIGNATED  BY  THE  DPERATOR.IT* 
REDUCES  THE  DATA  TO  CALCULATIONS  OF  NET  • 
POWER,  RAYLEIGH  AND  NUSSCLT  NUMBERS.  THE* 
UNCERTAINTY  ANALYSIS  IS  ALSO  INCLUDED.  • 


UARIABLES  USED  ARE: 

EMF  :  UDLTAGE  FROM  THE  THERMOCOUPLES. 

POWER  :  POWER  DISSIPATED  BY  THE  HEATERS 
T( I  ),  :  TEMPERATURE  CONUERTEO  FROM  THERMO¬ 
COUPLE  UOLTAGE 

Tavg  :  IS  THE  AVERAGE  TEMPERATURE  OF  THE 
CHIP.  IT  IS  obtained  MULTIPLYING 
THE  TEMPERATURE  FOUND  IN  EACH  FACE 
BY  THE  AREA  ANO  DIVIDING  BY  THE  TO¬ 
TAL  AREA 

Ts  :  CHIP  BACK  SURFACE  TEMPERATURE 
Tfiln  :  FILM  TEMPERATURE  OF  THE  DIELECTRIC 
Qnet  :  ELECTRIC  POWER  MINUS  CONDUCTION  LOSSES 
Tsink  :  AVERAGE  OF  THE  3  THERMOCOUPLES  IN 
THE  UPPER  HEAT  EXCHANGER 
Nul  :  LENGTH  BASED  NUSSELT  NUMBER 
Nu2  :  AREA-PERIMETER  BASED  NUSSELT  NUMBER 
D...  :  UNCERTAINTY  OF  A  VARIABLE  (EXCEPT 

D1 Iq  end  Celt) 

OTHER  VARIABLES  ARE  SELF-EXPLANATORY 


COM  /Co/  D(7) 

I 

DIM  Emf <  76  )  ,Pouer( 9  )  ,T( 76 ) ,Tavg( 9  )  .Tsi 9  ) 

DIM  Tf  il«i(  9  )  ,Onet<  9  )  ,H(  9  )  ,K<  9  )  ,Rho<  9  )  ,Cp(  9) 

DIM  N(9),Nu1(9).Ra1(9),Deit(9),Alfa(9),Pr(9) 

DIM  6r1(9>,Bete<9)  ,Dpouer( 9  )  ,Ra2( 9  ) 

DIM  6r2(9),Raf t(9),Raf2(9),Nu2(9) 

DIM  Roural ( 3  )  ,Rounu1 ( 3  ) 

I 

I  CORRELATION  FACTORS  TO  CONVERT  Enf  TO  DEGREES  CELSIUS.  SOURCE: 

I  HP  APPLICATION  NOTE  290,  P.  8,  NBS  POLYNOMIAL  COEFFICIENTS  FOR 
I  TYPE  T  (COPPER-CONSTANTAN )  THERMOCOUPLES. 

DATA  0.10086091 ,25727.9,-767345.8.78025596. 
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520  DATA  -9247485589.6.98511  .-2. 66E13. 5. 94E14 
530  I 

540  READ  D( •  ) 

550  I  PRECISION  RESISTOR  OALUE  IN  OHMS 
580  Rp-2.0 

570  1 

580  PRINTER  IS  701 
590  BEEP 
500  BEEP 
510  I 

620  INPUT  ‘ENTER  THE  NAME  OF  THE  FILE  CONTAINING  DATA"  .Ci df i leS 
530  I 

540  PRINT  USING  "10X.'‘''THE  RAW  Er,f  DATA  ARE  FROM  THE  FIlE:  “  1  0A '  :  0 1  df ;  i  e$ 

650  1 

660  INPUT  "ENTER  THE  °OUER  SETTING  ", Power? 

670  print  using  “9X."“  THE  POWER  SETTING  PER  CHIP  UAS:  " "  l0A"!Pouer» 

680  I 

690  INPUT  "ENTER  THE  TYPE  OF  LIQUID  USED "  ,L i qui d$ 

700  print  USING  "I0X.""THE  FLUORINERT  USED  wAS:  “ " , 1 0A" ; L iqu i dS 
710  input  "ENTER  THE  '7PE  OF  DIELECTRIC : 0=FC-75 , 1 -FC-43  ;=FC-7r’  Dlic 

720  I  .  ,  M 

730  INPUT  “ENTER  THE  WALL  SPACING" ,Wal IS 

740  print  USING  "10X."“THE  DISTANCE  TO  THE  FRONT  WALL  WAS:  ""  10A"!Wa"$ 

750  !  '  " 

760  INPUT  ‘ENTER  THE  GEOMETRY  TYPE:  0-HORIZONTAL  ,  1 -UERTICAL"  , Geo 
770  INPUT  ‘ENTER  THE  CHIP  ORIENTATION"  ,Chip$ 

780  print  USING  ''10X,""THE  CHIP  ORIENTATION  WAS:  ""  10A"!ChiD$ 

790  I 
800  BEEP 
810  BEEP 

820  ASSIGN  @File  TO  OldfileS 
830  ENTER  @FiieiEnf(»  ) 

840  I 

850  I . . . . 

860  I  CONVERT  Ernf  TO  DEGREES  CELSIUS  • 

870  . . . . . 

880  I 

890  FOR  >0  TO  60 

900  Sui*i-0 

910  FOR  J-0  TO  7 

920  Sun»Sum+0(  J )»Enf  (  I )" J 

930  NEXT  ; 

340  T < I  )“Sun 

950  NEXT  I 

960  FOR  1-71  TO  76 

970  Sun-0 

980  FOR  J-0  TO  7 

990  Sun-Sum+D( J  )»Enf < I  )“ J 

1000  NEXT  J 

1010  T(I)-Sun 

1020  NEXT  I 

1030  I 

1040  . . . 

1050  I  CONVERT  Enf  TO  POWER  * 

1060  I . . 

1070  I 
1080  J-1 
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t090  yolt-Enf(61 ) 

1100  FOR  I-G2  TO  70 

1110  Poyer<  J )-Enf(  I )• ( Uol t-Enf < I  )  )/Rp 
1120  J-J+I 
1130  NEXT  I 
1140  I 

I  I  50  Seo-0  THEN 
1160  L1-8.E-3 

1170  Alef-4.8E-5 

1 180  AriO”4.8E-5 

1190  Atop-1.44E-4 

1200  Abot-1.44E-4 

1210  Le-L1»1000 

1220  PRINT  USING  "  1  0X  "HORIZONTAL  LENGTH  SCALE  IS  (MM):  ,DD.D''iLe 

1230  ELSE 

1240  L1-2.4E-2 

1250  Alef-1.44E-4 

1 260  Arig« 1 . 44E-4 

1 270  Atop“4 . 9E-5 

1290  Al3ot"4.8E-5 

1 290  Le“L 1 • 1 000 

1300  PRINT  USING  1  0X  ,  “UERI'ICAL  LENGTH  SCALE  IS  (MM):  ■"•.DD.C":L: 

1310  END  IF 
1 320  Acen- 1 . S2E-4 
1330  Atot“5.7GE-4 
1540  t 

1350  I . . . * . 

1360  'CALCULATE  THE  AUERAGE  TEMPERATURES  OF  THE  BLOCK  FACES 
1370  I 

1 380  Tavg( 1  )«(T(0)»Acen+T( 1  )«Atop+T(2)»Ario+T(3)»Alef  +  T(4  >»Abct  l/Atot 
1390  TavQ(2)“(T(6 )»Acen+T (7 )*Atop+T(8 )«Arig+T (9)»Alef+T( 10)»Abot  >/Atot 

1400  Tavg(3 )“(T(  12 )»Acen+T(  13 )*Atop+T(  15  >»Arig+T( 16  l^Abot  +  Alef^Tl 14  )  )/Atot 
1410  Tavg(4)“(T(16)*Acen+T(19)*A'top+T(20  )*Arig+T  (21  )«Aief  +  T(22>»Abot  )/Atot 

1420  Tavg(5)“(T(24)*Acen+T(25)»Atop+T(26)»AriQ+T(27>»Alef+T(26)»Abo't)''Atot 

1430  Tavg(6  )-( T( 30  )*Acen+T(3l  )*Atop+T(32  )«Arig+T( 33  )»Alef  +  T( 34  )»Abot  )/( Atot  ) 
1440  Tavg( 7 )“(  T(  36 )»Acen+T  (  37 )»At op+T  (38)»Arig+T(39>»Alef*T>  40)«Abot  l/Atot 
1450  Tavg(6)“(T(42)«Acen+T(43)»Atop+T(44)*Arig+T(45>«Alef  +  T(46>»Abot  )/(Atot  ) 

1460  Tavg( 9  )“( T( 48  )*Acen+T( 49  )«At op+T( 50 )»Arig+T( 51  >»Alef f T( 52  )*Abot  l/Atot 
1470  I  ACCURACY  OF  THERMOCOUPLES  IS  0.5  DEGREES  C. 

1480  Dt-.5 

1490  Dtavg«.5 
1500  Dtsink“.5 
1510  I 

1520  I  RESISTANCE  AND  UNCERTAINTY  OF  PLEXIGLASS  FOUND  WITH 
1530  !  A  CONDUCTIVITY  OF  0.195  W/n.K  &  A  LENGTH  OF  12.0  MM 
1540  Rc-320.51 

1550  Drc-10.05 

1560  1 

1570  I . 

1580  1  CHIP  BACK  SURFACE  TEMPE^ ATURES  • 

1590  . . . . . 

1600  Ts( 1  )-T(5  > 

1610  T8(2)-T(11) 

1620  T9(3)-T(17) 

1630  T9(4)-T(23) 

1640  T8(5)-T(29) 

1650  T9(6)-T(35) 
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1660  Ts(7  )«T(41  ) 

1670  Ts(8)-T(47) 

1680  Ts(S)-T(B3) 

1690  Ts5uff“0 
1700  FOR  J-1  TO  9 
1710  Tssun^Tssurn  +  Ts (  J  ) 

1720  NEXT  J 
1730  I 

1740  Tsavg-Tssufi/S 
1750  I 

1760  I . . 

770  I  CONDUCTION  LOSS  CALCULATION  • 

1780  I . . . 

1790  Q1o5s1-(T'5)-T<57>)/Rc 

1800  01os52-(T(  1  1  l-TvbS  ) /.'Rc 

1810  01os53“( T< 1 7  )-T( 59  )  )/Rc 

1820  Oloss4>=(T(23)-TC60))/Rc 

1830  01os55=(T(29)-T(71  ))/Rc 

1840  Q1o5sE-(T(55  )-T( 72  )  )/Rc 

1850  01o557-(  T(  41  )-1  ( •  3  '  )/rxc 

1860  01oss3=(T(47)-T(74  .  )/R, 

1870  Oloss9*< T( 53  i-T( 75  !  )/Rc 

1880  01c55  =  (01o5El+Qloss2  +  QlossZ*0icss4fCloss5+01os6B+QIoE47-CIoss6'»0IcsES)''5. 

1890  DqloES'( Qlos:  )♦( (Dt /T(57  )'2  )  +  (Drc/Rc  I'C  : '  .5 
1 900  1 

1910 

1920  I  AUERAGE  SINK  TEM=ERATURE  CALCULATION  * 

1930 
1 940  I 

1950  TsinK-(T(54  )+T(55  .'  +  T(56  )  >/3. 

1 950  I 

i  970  PRINT  USING  "  1  0X  "AUERAGE  SINt-  TEMPERATURE  '■  C  / :  ''  ,DD2  .  D2  '  s  T  s  i  nt 

’  980  I 
'950  PRINT 
2000  I 

2010  I  TWO  CHARACTERISTIC  I.ENGTHS  WILL  BE  USED  TC  CAlCUl^TE  NUSSElT  NUMBERS; 
2020  I  Lt  IS  BASED  ON  THE  UERTICAL  OR  HORIZONTAL  DIMENSION  OP  THE  CHIP 
2030  I  L2  IS  BASED  On  THE  SUMMATION  OF  THE  AREAS  DIUIDE:  S';  THE  F  IMETER 
2040  I 

2050  L2-(2.»(6.»24./60.  )+2 . * ( 6 . *6 . /2e .  )+8.»24./B4.  )».00) 

2000  ! 

2070  PRINT  USING  ’Sv/'CHIP  Or.et(W)  Tavg-Ts  Nu’  Nu2  XUNC  IN  N 

U" "  ,10A” 

2080  PRINT 
2090  ! 

2100  . . . . . . . 

2110  I  CALCULATION  OF  NET  POWER.  Nu ,  Ra  AND  UNCERTAINTIES  • 

2120  I . . . 

2130  ! 

2140  FOR  J-1  TO  9 
2150  1 

2160  I  Dpouer  IS  BASED  ON  ACCURACY  0"  THE  VOLTMETER  AND  THE  PRECISION  RESISTORS 

2170  I  Dv  FOR  THE  VOLTAGE  DROPS  ACROSS  THE  CHIP  HEATERS  OR  PRECISION  RESISTORS 
2180  I  IS  5E-B  V.  Drp-0.05  OHMS. 

2190  DV-5.E-G 
2200  Drp-.05 
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2210  Dpouer(  J  l-PowerC  J  )•(  (Dv/Enf  (  J+El  )  )''2+(  Dv/(  L/ol  t-Enf  (  J+Bl  ) )  )''2  +  (Drp/Rp  )*2  )*  . 
S 

2220  ! 

2230  !  CALCULATION  OF  Qnet 

2240  Onet < J  )“Power( J  )-01o5S 
2250  Dqnet«(Dpower( J  )*2+Dq1oss*2  I' .5 
2260  ! 

2270  !  CALCULATION  OF  Tfiln 

2280  Tf i ln( J  )•( Tavgl J  l+Tsi nk  )/2 
2290  I 

2300  I  CALCULATION  OF  A  DELTA  TEMPERATURE 
2310  De 1 1 ( J  )“Tavg( J  )-Ts ink 
2320  Ddelt«(DtavQ*2+Dt6ink-2  )',5 
2330  I 

2540  I  CALCULATION  OF  CONUECTION  COEFFICIENT 
2350  H< J  )«Qnet ( J  )/( Atot •Del t ( J  )  > 

2360  Dh«H(  J  )♦(  (  Dqnet/Qnet  (  J  )  )'2+<DdeU/Delt  (  J  )  )“2  )"  .5 
2370  ! 

2380  i  PHYSICAL  PROPERTIES  ARE  TAKEN  FROM  THE  1965  3M  PRODUCT  MANUAL 
2390  !  FOR  FLUORINERT  ELECTRONIC  LIQUIDS 

2400  I 

2410  IF  DUq-0  THEN 
2420  ! 

2430  I  calculation  OF  FC-75  THERMAL  CONDUCTIUITY 
2440  K( J  )-( .65-7.89474£-4»Tfiln( J  ))/10 
2450  ! 

2460  I  CALCULATION  OF  FC-75  DENSITY 
2470  Rho(  J )-{  1 .825-. 00246 •Tfiin( J  ))*100e 
2480  1 

2490  !  calculation  OF  FC-7S  SPECIFIC  HEAT 

2500  Cp( J  )■(  .241 n l+3.7037E-4«Tf iln( J ) )»4187 
2610  I  THE  4187  CONUERTS  FROM  CALORIES  TO  JOULES 
2520  I 

2630  I  calculation  FC-75  KINEMATIC  UISCOSITY 

2540  N<  J  )=1 .4074-2. 964E-2»Tfi In' J  )+3 . 801 8E-4»T f i In' J )'2-2 . T30SE-6«Tf i ln( J  i' 3+8. 
1679E-9»Tf  ;in(  J  )''4 

2550  I  CONVERT  FROM  CENTI5T0KES  TO  n"2/s 
2560  N( J  )-N( J  )*1 .E-E 
2570  I 

2580  !  CALCULATION  0-  THE  COEFFICIENT  OF  THERMAL 

2590  I  EXPANSION  IBETAI 

2600  Bet  a;  J  )•  .  00246  M  .  62E-  .00246  *Tf  i  In'  J  )  ) 

2610  I 
2620  END  IF 
2630  I 

2640  IF  Dliq«I  THEN 
2650  I 

2660  I  calculation  OF  FC-43  THERMAL  CONDUCTIVITY 
2670  K( J  >-(  . 0666-9. 864E-6«Tfiln(J  )) 

2680  I 

2690  I  CALCULATION  OF  FC-43  DENSITY 
2700  Rho(J  I .913-.00218*Tfiln( J))»1000 
27)0  I 

2720  I  calculation  OF  FC-43  SPECIFIC  HEAT 
2730  Cp( J  )-( .2411 1 1+3.7037E-4«Tfiln( J))*4167 
2740  I 

2750  I  calculation  OF  FC-43  KINEMATIC  VISCOSITY 
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2760  N(  J  )-8. 875-.  47007*7  filn(J  )+l  .387E-2»Tf  iin(  J  )'2-2  .  UB9E-4*Tf  i  1p-(  J  )"3+1  .3139 
E-B»Tfiln( J  )'4 
2770  I 

2780  N< J )-N( J  )•! .E-6 
2790  ! 

2800  I  CALCULATION  OF  COEFFICIENT  OF  THERMAL 
2810  I  EXPANSION  [BETA] 

2820  Beta( J  )«.00218/(1 .913-.00218*Tfiln( J  ) ) 

2830  I 
2840  END  IF 
2850  I 

2860  IF  DliQ-2  THEN 
2870  ! 

2880  I  CALCULATION  OF  FC-71  THERMAL  CONDUCTIVITY 
2890  K(J)-. 71/10 

2900  I 

2910  I  CALCULATION  OF  FC-71  DENSITY 
2920  Rho(  J  )-(2.002-.00224«Tfili«i(  J))«1000 
2930  I 

2940  I  CALCULATION  OF  FC-71  SPECIFIC  HEAT 
2950  Cp< J  )“(  .24  11 1 l+3.7037E-4*TfiIn( J))*4187 
2960  I 

2970  I  CALCULATION  OF  FC-71  KINEMATIC  VISCOSITY 

2980  N( J )-EXP<6.8976-. 1388»Tfilm( J  )+l .331E-3*Tfilm( J  r 2-7 . 04 1 E-6*Tf 1 1 J  ) '3+ 1 . 5 
23E-8»Tf 1  In! J  )'4  ) 

2990  ! 

3000  N(  J )-N< J  )*1  .E-B 
3010  I 

3020  I  CALCULATION  OF  THE  COEFFICIENT  OF  THERMAL 
3030  I  EXPANSION  [BETA] 

3040  Beta! J  )». 00224/ ( 2 . 002- .00224»Tf i ln( J  )  ) 

3050  ! 

3060  END  IF 
3070  I 

3080  '  calculation  OF  THERMAL  DIFFUSIVITY  [ALPHA] 

3090  Alfa! J  )«K! J  )/<Rho!J  )»Cp< J  )) 

3100  I 

3110  1  CALCULATION  OF  PRANDTL  NUMBER 

3120  Pr( J  )=N( J  )/Ai fa! J  ) 

3130  I 

3140  I  CALCULATION  OF  NUSSELT  NUMBERS 
3150  Nul ( J  )«H( J  )*L1 /K( J  1 
3160  Nl2!  J )-H( J  )»L2/K! J  ) 

3170  Dnul-NuK  J  )*!Dh/H(J  )) 

3180  Pernu 1 “( Dnu 1 /Nu 1 lJ  )  )• 1 00 
3190  I 

3200  '  calculation  OF  6RASH0F  NUMBERS 

3210  6rt  (  J  )-9.81  ‘Beta!  J  )•!  LI  '3  )*Delt  ( J  )/N(  J  )'-2 
3220  Gr2' J  )  =  5.ei ‘Beta! J  )•! L2'3  )*Del t( J  >/N( J  r2 
3230  Dc"  1 “G' I ( 2  ^ ( Doe i t /Del 1 1 J  ^ ; 

3240  I 

3250  I  CA.CULATION  of  RAYLEIGH  NUMBERS 
3260  Re  1 ! J  l“Gr H J >*Pr< J  )♦ 1 . E-6 
3270  Ra2( J  )-Gr2i J  )»Pr( J  )•! .E-E 
3280  D’"e  1  "Re  1  (  J  )  •  (  Doi"  1 /Sr  1  (  J  )  ) 

3290  Perra I “( Dra 1 /Ra 1 ! J  1  )• 1 00 
3300  I 
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3310  I  CALCULATION  OP  FLUX  BASED  RAYLEIGH  NUMBERS 

3320  Raf 1 <  J  )-( ( 9.81 ‘Beta; J  )*L1 '4»Qnet( J ) )/(K( J  )»N( J  >»Alfa(  J )*Atot ) >•! .E-B 
3330  I 

3340  Raf2( J  )-( ( 9.81 ‘Beta <  J )«L2*4*0net(  J ) >/(K(  J )*N(  J '^Al f a'  J )»Atot ) )* 1 .E-E 
3350  I 

33B0  I  DATA  AND  UNCERTAINTY  OUTPUT 
3370  I 

3380  PRINT  USING 
2 ( J  )  .Pernu 1 
3390  PRINT 
3400  I 

3410  PRINT  USING 
<  J  > 

3420  PRINT  USING 
1(J) 

3430  PRINT  USING 
3440  PRINT 
3450  PRINT  USING 
DO. 00" iOnut 
34G0  PRINT  USING 
D.DDD" iRa1  (  J ) 

3470  PRINT  USING 
DD.DD" iDral 
3480  PRINT  USING 
DD.DD” (Parra  I 
3490  PRINT 
3500  NEXT  J 
3510  I 

3520  Ralsun“0. 

3530  Raf2sun«0. 

3540  Nul5un“0. 

3550  Nu2sun"0. 

35B0  Onetsun“0. 

3570  Delt5un«0. 

3580  FOR  J-1  TO  9 
3590  Ra  1  SLin"Ra  1  <  J  )+Ra  1  sun 
3G00  Raf2sun“Raf2(  J  )-*-Raf2sun 

3B  1  0  Nu  1  SLin“Nu  1  (  J  )+Nu  1  sun 
3B20  Nu25un“Nu2 ( J  )+Nu25un 

3B30  Onet 8un“0net ( J  l+Onetsun 

3B40  Delt sun“De  1 1  <  J  )  i-De  1 1 sun 

3B50  NEXT  J 

3BB0  Ra 1 ( Avq  )“Ra 1  sun/9 . 

3E70  Raf2< Avq  )“Raf2sun/9 . 

3B80  Nu 1 ( Avq  )“Nu I  sum/9 . 

3G90  Nu2( Avq  )"Nu28un/9 . 

3700  Onet ( Avq  )“One t sun/ 9 . 

3710  Delt( Avq >“Deitsun/9. 

3720  I 

3730  FOR  J-!  TO  3 
3740  Roura I <  J  )'0 
3750  Rouinul<J)“0 
3760  NEXT  J 
3770  FOR  J-1  TO  3 

3780  Row-al ( J )-(Ral ( J  >+Ra1 ( J+3  )*Ral ' J+B  )  )/3.0 
3790  Rounul(J)“(Nul(J)  +  Nul(J  +  3>'*'Nul(J-*-B))/3.0 
3800  NEXT  J 


“  10X .0 .IX .4(  4X .DDD.DD . ) ,6X .ODD. DD" iJ  .Onet ( J  )  .Delt( J  )  .Nul  ( J  )  ,Nu 

■12X."”TEMP  BASED  RAYLEIGH  NUMBER  *  E-E  IS:  .DDDDDD.DDD" (Ral 
■12X.'"FLUX  BASED  RAYLEIGH  NUMBER  •  E-B  IS:  *" .DDDDDD.DDD" (Raf 
•  1  2X  . "  AUERAGE  TEMPERATURE :  “  ‘  .ODD .  DD"  :  Ta vqI  J  > 

•BX.-“UNC  IN  THE  NUSSELT  NUMBER  (Nul)  IS: 

"EX. ""TEMP  BASED  RAYLEIGH  NUMBER  •  E-B  IS:  *“.DDD 

“BX,""UNC  IN  THE  TEMP  BASED  RAYLEIGH  NUMBER  •  E-B  IS: 
"GX.""5iUNC  IN  THE  TEMPERATURE  BASED  RAYlEIGH  NUMBER  IS:  "".D 
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3810  I 
3820  PRINT 
3830  PRINT 
3840  PRINT 

3850  PRINT  USING  ■' 1  ZX  "ARRAY  AUG  Ral*E-6  IS  : "  “  .DDDDD .  ODD"  ;  Ra  H  Avg  ) 

38B0  PRINT  USING  " 1 2X  "ARRAY  AUG  Nul  IS : ‘ "  ,DOD . DD’ sNu 1 ( Avg  ) 

3870  PRINT 

3880  PRINT  USING  "12X."'ROU  1  AUG  Rat«E-6  IS: .DDDDD.DDD’ :Roural ( 1) 

3890  PRINT  USING  "I2X,""ROW  1  AUG  Nul  IS: “ " .DDD.DD" ;RownuH  1) 

3900  PRINT 

3910  PRINT  USING  "l2X.""R0U  2  AUG  Ral*E-B  IS : ‘ ' .DDDDD . ODD" ( Roura 1 ( 2  ) 

3920  PRINT  USING  "12X,"“R0U  2  AUG  Nul  IS: " ‘ .DDD .DD" i RownuH 2  ) 

3930  PRINT 

3940  PRINT  USING  “12X,"“R01J  3  AUG  Ral«E-6  IS:  .DDDDD.DDD"  :RouraH  3  ) 

3950  PRINT  USING  "12X.“"R0U  3  AUG  Nul  IS: " " .DDD.DD" ;Roynul( 3  ) 

39G0  PRINT 

3970  PRINT  USING  " 1 2X  "ARRAY  AUG  Qnet  IS:  ""  .D .DD" i Qnet ( Avg  ) 

3980  PRINT  USING  " 1 2X  ARRAY  AUG  <TAU6-TSINK)  IS:  “ "  .DD . DD" :Del t ( Avg > 

3990  ASSIGN  CFile  TO  • 

4000  END 
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APPENDIX  C.  COMPtTTER  PROGRAM  ACQ2 


ic 

^  RI_E  ACC; 

:e 

( 

5C 

'  uiR:~'£k  E'  ^CDR 

4e 

se 

THIS  program  AO 

Be 

1  FROM;  THE  HP3456 

7e 

'  POP  THE  3  E'  3 

80 

'  IT  IS  URITTEk  I 

se 

DUM  -  DIGITA^ 

:0e 

'  DAS  -  data  ACOU 

I  !C 

120 

DIM  Ei'^  *  64  )  .Pouie- 

i3e 

1 

1  40 

'PRECISION  RESIST 

150 

Rp*2 . 0 

ise 

1 

170 

1  ASSIGN  OUTPUT  T 

180 

PRINTER  IS  701 

1  90 

BEER 

200 

input  >£nteR  THE 

210 

1 

220 

IF  In»t  then 

230 

BEEP 

240 

INPUT  -enter  THE 

250 

260 

PRINT  USING  ■■  15X  , 

270 

else 

260 

BEEP 

2  90 

input  "ENTER  THE 

300 

PRINT  USING  "lOX, 

310 

END  IF 

320 

PRINT 

330 

PRINT  USING  -15X ^ 

340 

1 

350 

INPUT  "ENTER  the 

360 

PRINT  USING  "IBX 

370 

1 

380 

input  "ENTER  the 

350 

PRINT  USING  "IBX 

400 

1 

410 

input  'ENTER  THE 

420 

PRINT  USING  "15X 

430 

IP  iri-i  THEN  ass: 

440 

1 

450 

IP  Ii«e0  THEN 

460 

CREATE  BOA*  NewT 

470 

assign  PPile  TC  1 

480 

End  ir 

490 

500 

1  READ  DATA  INTO 

:  'H:i*1=SDN  5- • 

litres  the  UOl’AdE 
f  DL'K  L'lft  the  HPZtS~!h  DAS 

,4,CRAV  NU'SC  CIRCUIT  BOARD, 
k  BASIC  :.e. 

'Catheter 


-.Old' 


.  1  PA'  i  Nci^f  ; 


,  10A" ;Wol IS 


I  READ  DATA  INTO  H09E2B  COMPUTER.  70S  IS  THE  DAS.  122  IS  THE  DUM. 


lies 


ieS 
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5ie  I 

520  IF  Ir"»e  THEN 

530  I  AP  RESETS  DAS.  AF  IS  FIRST  CHANNlEi- .  AL  IS  LAST  CHANNE_. 

540  0l''°UT  705;  "AR  A^OO  AlIZ-" 

550  '  F'  SETS  FUNCTION  TO  DC  UOL'^S .  RI  SETS  RANGE  TO  AU^O. 

560  I  T1  SETS  TRIGGER  TO  INTERNA^.  20  SETS  AUTO  ZERO  TO  OFF. 

570  I  FLO  SETS  FILTER  TO  OFF. 

590  OUTPUT  7::;  "FI  R1  T 20 
ESC 

600  '^OR  1=0  TO  :Z 

E10  I  AS  CAUSES  THE  DAS  TO  ANALOG  STER  THROJG^^  THE  ChANNR..= 

620  OUTPUT  705; "AS" 

E30  WAIT  1 

640  I  ENTER  SENDS  UO_’AGES  FPOr  DOT  TO  DAS. 

ESC  ENTER  722  :Er/i  I 

660  BEER 

670  NEXT  I 

680  OUTPUT  70S;"AF  A^R ’  ALE4  " 

ESC  POP  1-41  TO  64 
700  OUTPUT  705: “AS' 

7:0  WAIT  1 

720  ENTEF  722;Erfv I  ; 

730  BEEP 

740  next  I 

750  OUTPUT  @File;Er.f  (  • 

750  I 

770  Else 

780  ENTER  eFile:En*C  ) 

790  END  IF 
B00  I 

810  I  AR  RESETS  DAS 
620  OUTPUT  70St"AR" 

850  I 

840  I  CONUERT  THERMOCOUPLE  U0LTA6ES  TO  TEMPERATURE  IN  DEGREES  CELSIUS. 

850  I  THERMOCOURLES  CALIBRATED  AGAINST  PLATINUM  RESISTANCE  THERMOMETER 

880  I  MARCH  2-3,  1S82.  TEMP.  RANGE  10-100  DEG.  C. 

970  t  calibration  CURUE  fit  BY  TABLECURUE  SOFTWARE. 

880  I  NOTE:  CAL.  CURUE  UOLTAGES  ARE  MILLIUOLTS,  Em'lI'S  ARE  IN  MICROUD"' 

890  A=. 24577463 

900  8-24.898086 

910  0=-. 075219169 

520  FOR  I- 11  TO  13 

530  T ( I  )«A+B«Enf ( I  )♦ 1 .E*3+C»( 1 .E+3»Enf ( I  ) )'  3 

940  NEXT  I 

950  FOR  1=41  TO  59 

360  T(  I  )«A+B»Er>)f  <  I  )*1  .E*3+C*(  1  .E*5»Emf(  I  )  )'3 
970  NEXT  I 
980  I 

990  I  CHIP  TSE'S  (THE  TRANSISTORS)  ALSO  CALIBRATED  SAME  TIME  AS  ABOUE. 
1000  I  TEMP.  (DEG.  C>  US.  BASE-EMITTER  (Ube)  UOuTAGE  CALIBRATION  CURUER 
1010  I  PRODUCED  USING  TABLECURUE  SOFTWARE. 

1020  T(60  )-577. 58074-575. 54353»Er.f(60) 

1050  T(61  ’=575. 86885-574. 05057*Enf(Gl ) 

1031  T(6:  )-577. 64556-575. 50862*Emf(E2 ) 

1022  ..c-’E, 96331-574. £S738»Et',f'.E3) 

1033  T(64  )“57e. 29580-575. 9e7S3*Emf(64 ) 

1050  ' 


88 


1  050 

pp:^■ 

I'EINE 

’ 5>  , ’ " VOl 

TRGE  SJPPlV  WHS:  “ '  ,02  .  D2 '  :  Ef"' '  1  0 

1070 

PR  IK" 

UEIN5  '■ 

^5>'  ‘  ftME 

lENT  TEMP  WHS:  " "  .DC .  D '■  i  T  50' 

1080 

FFIN" 

1090 

'  POiJ 

E  ^  C  ^  -  C  U 

.f’ZOr.S 

1  1  00 

i 

•  1  10 

“  0^  I 

=  ’  TC  3 

■  ';c 

=  :  -r 

r  .  *  r  ^  z 

"♦P; 

-42 

k5  ■  " 

J  Si. 

■ 

■  ’50 

r  "lE  * 

=  '■’■0  3 

1  1B0 

■  2  ’^E * 

'  1  >  ♦  E  ^ '  r 

-"5  >  >9t 

!  •  70 

1  1  62 

NE> 

I 

1  150 

=^0F  I 

=7  TC  5 

1Z00 

Roue 

- : I  >=Er‘ 

'.4  >  ♦  E  f  ^  5 

1 '  ■■  3 ‘Re  ' 

':’0 

KE>'" 

. 

i::0 

’ 

1130 

SEER 

1340 

' 

1250 

PR  IN'" 

l'EIN'3  “ 

ie.>  ,  ■ '  R__ 

tehpERH"UP£S  HPE  in  DEGREES  CELSIUS . 

1250 

1 

1270 

“PIN" 

1280 

1 

12  90 

ERIN" 

U: I N5  ' 

1  2>  ,  ■  2-‘I 

P  CP'IP  Li:  POwERi  U- 1" '  ' 

1300 

PRINT 

1310 

PRIN7 

1220 

PPIN7 

USING  " 

IX  /'"CRIP 

NOl:  N/H  N'H 

,D.DDD" ;Powe 

r(  1  ) 

1230 

PRINT 

USING  ■■ 

12X 

( HUG ; . 

1340 

PPIN-^ 

1350 

PRINT 

USING  - 

IX  /  "CHIP 

N02:  "•  .2(DDD.DD.10X  ;  ,C.DDC"!T(B0 ' ,T(45 ) .Power (2 ) 

1360 

PRINT 

USING  •' 

1  2X  /  " 

( HUG  ) . 

1  370 

PRINT 

1380 

PRINT 

USING  " 

IX," "CHIP 

N03:  N/H  N/H 

.0 . DDD ' i Powe 

1390 

PRINT 

USING  •' 

12X  " 

(  HUG  )  " 

1  400 

PRINT 

1410 

PRINT 

USING  " 

IX," "CHIP 

NC4:  "  ••  .2'DDC.DC  ,iex  '  .C.DDC"  :7>-  El  1  .T(  4E 

' .Powe- ' 4  ' 

1420 

PRINT 

1430 

PRINT 

1  440 

PRINT 

USING 

IX  ,”''CHIP 

NOE:  •' '  .2(DDD.DD  ,  10X.  '  ,D.DDD'  ;T(  52  ,T(  47 

.Power  (  5  > 

1450 

PRINT 

1460 

PRINT 

1470 

PRINT 

USING  ■■ 

IX," "CHIP 

NOE:  .2' DDD.DD  ,  1  0X  )  .D  .  ODD" ;  T  (  53  I  .T  (  46 

1  .Powe- ( 5 ) 

1480 

PRINT 

14  90 

PRINT 

1500 

PRINT 

USING  " 

IX "CHIP 

N07:  N/H  N/H 

.0 . DDD" ; Powe 

r(  7  ) 

1510 

PRINT 

U'SING 

1  2X  , " " 

( HUG  )" " " 

1520 

PRINT 

1530 

PRINT 

USING 

IX .""CHIP 

N06:  •"  .2(DDD.DD  .10X  )  .0 . DDD" ; T( 84 ) .T : 49  > .Powe-i  8 ' 

1540 

PRINT 

USING  ” 

12X 

( HUG  >" " ' 

1550 

PRINT 

1560 

PRINT 

USING  ” 

IX." "CHIP 

N09:  N/H  N/H 

.0 . DDD" :Powe 

r'  9  } 

1570 

PRINT 

USING 

12X  ,"  " 

(HUG  )■  "• 

1580 

PRINT 

89 


tE!^=E~^’UP£S 


’E9C 

=-IN’  L'EINC- 

■'  1  2 

"CIRCUIT  E 

lOARD  SUESTRATE  SURpaOE  tem-EPatureS  " " ' 

isoe 

PRIM 

1£!0 

PPIN”  l'£!N0 

■’ !  y  ,  ■■  '■ 

p:  ""  ,DDC 

DD ■  ;  T ; 4  1  ' 

1E20 

PPIMT  L'SINC- 

•  1  y  , '  ■' 

S:  "'.DDD 

DO"  ;T; 42  ) 

1£3e 

UEZNS 

C:  "".DDD 

DD" ;T ( 43  ) 

1  £40 

PRIM'  USING 

2:  "".DDD 

DD" :Tk  44  : 

i£5e 

PR  IN’ 

1  5£0 

PRIN'  USING 

"  5  >  ,  ■  ' 

HEAT  exchanger  TEMPERATURES:  RIG-:'  CENTER 

1£70 

PR  IN'  USING 

"  1 0y . ' 

"BOTTOM  IE 

INSULATED . 

lEee 

PRINT  USING 

"  1 0  y , 

"TOP  :  "  ■■  .24 

tX  .3(  DD.DD  ,EX  '  i T-,  ’  ;  ,  ,T.  :  2  ■  ,t.  ■  2  i 

itse 

PR  IN' 

1700 

1 

1710 

PRINT  USING 

"SX  ' 

CIRCUIT  BOARD  ASSEMBLY  BACK  TEMPERATURES  ARE: . 

1720 

PRINT 

1730 

PRINT  USING 

1 0>: .  ‘ 

"CHIP  NOl 

"  .2X  .OD.DD"  ;T(  El  ) 

1740 

PRINT  USING 

"  1  0X  .  •' 

"CHIP  no: 

'•  .2X  .DD-DG" iT<52 ) 

1750 

PRINT  USING 

'  1 0x  .  ■■ 

'CHIP  N03 

■  •  ,2>:  .DD.DD"  ;T(S2 

17B0 

PRINT  USING 

'  1 0y. .  '■ 

"CHIP  N04 

■ "  .2X  .DD.DD" iT(54 ; 

1770 

PRINT  USING 

"10) 

"CHIP  NOE 

,2X  .DD.DD"  ;T!S£:- 

!7e0 

PRINT  USING 

'  1 0>:  .  '■ 

"CHIP  NOG 

*"  ,2X  .DD.DD"  :T(5£  : 

1790 

PRINT  USING 

" 1 0y  .  ’ 

"CHIP  NOT 

""  .2X  .DD.DD" ;T)57  ) 

1800 

PRINT  USING 

•'  1  C>  . " 

“CHIP  NOG 

.2X  .DD.DD"  iT(58  > 

1810 

PRINT  USING 

■■  1  C>  , ' 

"CHIP  NOS 

" "  ,2X  .DD.DD" iT(5E^ 

1820 

BEEP 

1850 

1  REASSIGN  PRINTER 

'0  the  Of 

■'T 

1  840 

PRINTEP  IS  1 

1850 

ftSSIGN  efile 

TO  ♦ 

1850 

EN2 
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APPENDIX  D.  COMPUTER  PROGRAM  CALC2 


10 
20 
30 
40 
50 
50 
70 
80 
90 
100 
1  10 
120 
1  30 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
t 

500 


’  file  CALC2 

I 

I  URITTENi  BY  LCDR  R.  THOMPSON  MAR  92. 

I  MODIFIED  early  APR  92.  BASED  ON  THE  PROGRAM 

'  calcdiel  written  and  modified  by  PAMUK  , 

I  benedict,  TORRES,  AYTAR ,  MATTHEWS  AND  THOMPSON. 

I 

I  THIS  PROGRAM  ANALYZES  THE  DATA  READ  FROM 
I  A  DESIGNATED  "ACQZ"  DATA  FIlE.  ACQ2  DATA  FIlES 
i  ARE  FOR  THE  3  BY  3  ARRAY  NWSC  CIRCUIT  BOARD. 

'  IT  REDUCES  THE  DATA  TO  CALCULATIONS  OF  NET  POWER,  RAYuEIGH 
I  AND  NUSSElT  number.  THE  UNCERTAINTY  ANAlVSIS  IS  ALSO  INCLUDED. 

!  IT  IS  WRITTEN  IN  BASIC  2.0. 

I 

I  UARIABlES  USED  ARE: 

I  EMF  :  DATA  FILE  UOLTA6ES. 

!  POWER  :  POWER  DISSIPATED  BY  THE  CHIP  RESISTORS  (W) 

I  T(I)  :  TEMPERATURES  CONUERTED  FROM  THERMO- 
!  COUPLE  UOlTAGES  (DEG,  C) 

I  Ts  ,  :  CIRCUIT  BOARD  ASSEMBLY  BACK  TEMPERATURE  (DEG.  C) 

I  Tfilrt  :  DIELECTRIC  FILM  TEMPERATURE  (DEG.  C) 

I  Tavo  :  AUG.  TEMP.  OF  THE  5  CHIP  TSE’s  (DEG.  C) 

I  Qnet  :  ELECTRIC  POWER  MINUS  CONDUCTION  LOSSES  (W) 

!  Tsink  ;  AVERAGE  OF  THE  3  THERMOCOUPLES  IN 
!  THE  UPPER  HEAT  EXCHANGER  (DEG.  C) 

I  Nu  :  LENGTH  BASED  NUSSELT  NUMBER 
I  Rc  :  THERMAL  RESISTANCE  FOR  CONDUCTION  (DEG.  P/W) 

ID...  :  uncertainty  OF  A  VARIABLE  (EXCEPT 

I  D1 iQ  AND  Delt  > 

I  OTHER  VARIABLES  NOT  SELF-EXPLANATORY  ARE  DEFINED  IN  THE  PROGRAM 

I 

DIM  Efnf(G4  )  . Power (  1  :  9  )  ,T(  64  )  ,Ts(  1  :  9  )  ,T5ide(  1  :  9  >  .T]  idi  1  :  9  ) 

DIM  Qnet  (  1  :  5  )  ,Del  t  (  1  :  9  )  ,H(  1  :  9  )  .Nu(  1  :  9  )  ,Ra(  1  :  9  ) 

DIM  Gr(l:9),01o55(l:9),Dh(l:9  /  ,Dpouer( 1:9) 

DIM  Rowgr( 1:3)  ,Rowra( 1:3)  ,Rounu( 1:3) 

I 

I  PRECISION  RESISTOR  VALUE  IN  OHMS 
Rp-2.0 

I  ASSIGN  OUTPUT  TO  HP  THINKJET  PRINTER 

PRINTER  IS  701 

BEEP 

I 

PRINT  USING  "10X,“OATA  TAKEN  BY  THOMPSON""’ 

INPUT  "ENTER  THE  NAME  OF  THE  FILE  CONTAINING  DATA’ ,01df i le$ 

I 

PRINT  USING  "10X,""THE  RAW  Enf  DATA  ARE  FROM  THE  FILE:  " " , 1 0A " i 01 df i le$ 

I 

INPUT  "ENTER  THE  APPROX.  POWER  SETTING  "  .Powers 

PRINT  USING  "9X,""  THE  APPROX.  POWER  SETTING  PER  CHIP  WAS:  10A’' (Power 
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510  INPUT  'ENTER  THE  TYPE  OF  LIQUID  USED"  .Liquids 

520  PRINT  USING  "t0X.""THE  FLUORINERT  USED  WAS:  " "  ,  I0A“ -.Liquids 

530  INPUT  'ENTER  THE  TYPE  OF  DIELECTRIC : 0-FC-75 , 1 -FC-43  ,2-FC-7 1"  ,D1 iq 

540  I 

550  INPUT  "ENTER  THE  WALL  SPACING"  .Ual IS 

560  PRINT  USING  "tOX/'-THE  DISTANCE  TO  THE  FRONT  WALL  WAS:  " "  .  I 0A " i Wal 1 S 
570  I 
580  BEEP 

590  ASSIGN  ePile  TO  OldfiieS 
600  ENTER  @File;Enf ( •  ) 

610  I 

620  I  CONVERT  thermocouple  VOLTAGES  TO  TEMPERATURE  IN  DEGREES  CELSIUS. 

630  I  THERMOCOUPLES  CA_IBRATED  AGAINST  PLATINUM  RESISTANCE  THERMOMETER 

640  MARCH  2-3,  tSS;.  TEMP .  RANGE  10-100  DEG.  C. 

650  '  CALIBRATION  CURVE  “^IT  BV  TABLECURVE  SOFTWARE. 

660  ‘  NOTE;  CAl.  CURVE  VOLTAGES  ARE  MILLIVOLTS,  EnfC.'S  ARE  IN  MICROVOlTS. 

670  A*. 24977463 

680  6=24.896088 

690  C--. 07921 91 ES 

700  FOR  >11  TO  13 

710  T(  I  )-A+B»Emf:  :  )•:  .E■^3■^■C•(  1  .E•^3•Emf(I  ))"3 

720  NEXT  I 

730  FOR  >4  1  TO  5S 

740  T(  I  )-AtB»Enf ;  I  )«  '  .E-f3+C»'  1  .E+3*Emf  (  Ill'S 
750  NEXT  I 
760  I 

770  I  CHIP  TSE’s  (THE  TRANSISTORS  1  ALSO  CALIBRATED  SAME  TIME  AS  ABOVE. 

780  •  TEMP.  (DEG.  Cl  VS.  BASE-EMITTER  ( Vbe  1  VOLTAGE  CALIBRATION  CURVES 

790  I  PRODUCED  USING  TABLECURVE  SOFTWARE. 

800  T(60  1-577. 58074-575. 54353«EMf( 60) 

810  T(61  1-575. 66889-574. 05057«Enif(61  ) 

820  T(62  1-577. 64556-575. 50862«EHf<62  1 

830  T(63  1-576. 98331 -574. 69738«Erif(  631 

840  T(64  1-578. 28560-575. 99793«Ernf(64  1 

850  Tavg-(T(60>T(61  l-t-T(  62  >T<  63  l-t-T<  64  1  1/5 

860  ! 

870  I  POWER  CALCULATIONS 
880  FOR  >1  70  3 

890  Pouer(I  l-Enf{01»Enf<51/(3»Rpl 

900  NEXT  I 
910  J-4 

920  FOR  >1  TO  3 

930  Pouer< J  l-Enf ( I  l»Enf < 1+5  1/Rp 

940  J-J+1 

950  NEXT  I 

960  FOR  >7  TO  9 

970  Pouer(I  l-Emf<41«Enf(9)/(3»Rpl 

980  NEXT  I 
990  1 

1000  I  BELOW  LENGTHS  AND  AREAS  BASED  ON  CHIP  MODELED  AS  A  SQUARE  WAFER. 

1010  I  ALL  DIMENSIONS  IN  h  OR  SQUARE  i*i. 

1020  L1-8.89E-3 

1030  Alef-1.65E-5 

1040  Ang-Alef 
1050  Atop-Alef 
1060  Abot-Alef 
1070  Alid-7.90E-5 
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1080  Le-L1»1000 

1090  PRINT  USING  "10X  .““LENGTH  SCALE  IS  ( nm  ) :  ”“.D.DD“:Le 
1100  I  Atot  IS  TOTAL  AREA  FOR  CONUECTION 
1110  At o t “Al e f +  Ar  ig  + At op  +  Abot +A i 1 d 
1120  I 

1130  I  CHIP  SIDE  TEMPERATURES.  ASSUMES  IT  IS  THE  AUERAGE  OF  THE  TSE  TEMP. 

11 40  I  AND  THE  LID  TEMP..  UITH  CHIP  TEMPS.  1.4.7  ASSUMED  TO  BE  EQUAL  AND 
1150  I  OHIP  TEMPS.  3.E.9  ASSUMED  TO  BE  EQUAL. 

1  1E0  I 

1170  Tsidev 1  )-( T( 46  )+T( 61  > >/2 
1180  Tsidei 2  )“( T( 45  )+T(G0  )  )/2 
1190  Tsidet 3  )«( T( 46  )  +  T( 63  )  )/2 
1200  T5ide( 4  )-( T< 46  )  +  T( 61  )  )/2 
1210  T£ide(5  )«(T< 47  )+T{ 62  )  )/2 
1220  Tside(E  )-=(T(48  )  +  T'63  )  )/2 
1250  Tsi de( 7  )-( T< 46  )+T( 6  I  )  )/2 
1240  Tside(  8  )•=(  T(  49  )+T(  64  )  )/2 
1250  Tside{9.'=(T(4S)  +  T(E3))/2 

1260  ' 

1270  i  CHIP  LID  TEMPERATURES.  SAME  ASSUMPTION  AS  CHIP  SIDES. 

1280  Tlid(  1  )  =  T(4G  ) 

1290  Tlid(2  )  =  T(4E) 

1 300  Tiidi3)  =  T(4c; 

1310  T1 id( 4  )=T( 4£ J 
1320  Tiid(5)“T(47) 

1330  Tlid(6  )=T(48  ) 

1340  Tiid(7)-T(46  ) 

1350  Tlid(8)-T(49) 

1360  Tlid( 9  )-T( 48  ) 

1370  I 

1380  I  CIRCUIT  BOARD  ASSEMBLY  BACK  TEMPERATURES 
1390  T5Sun-0 
1400  FOP  I-l  TO  9 
1410  T5( I  )“T( 1+50  ) 

1420  T9Sun“Ts5un+Ts(  I  ) 

1430  NEXT  I 
1 440  I 

1450  Tsavo“T55un/ 9 
1460  I 

1470  I  CONDUCTION  LOSS  CALCULATION 

1480  I  ONE  DIMENSIONAL  CONDUCTION  ASSUMED  THROUGH  THE  FOLLOWING: 

1490  I  4.0GE-4  m  (2. 016’)  OF  SILICON  (CHIP) 

1500  !  7.B2E-4  M  <0.030”)  OF  ALUMINA  (CIRCUIT  BOARD) 

1510  I  1.27E-4  tn  (0.005*)  OF  SILICONE  RUBBER 

1520  I  I.18E-2  n  (0.465”)  OF  ACRYLIC  (PLEXIGLAS) 

1530  I  UALUES  FOR  THERMAL  CONDMCTIUITY  K  IN  W/n-DEG.  C.  SOURCE: 

1540  I  SILICON.  ALUMINA:  NWSU .  MR.  TONY  BUECHLER 

1550  I  ALL  OTHERS:  MATERIALS  ENGINEERING  1978  MATERIALS  SELECTOR.  UOL .  86. 
1560  I  NO.  6.  REINHOLD  PUBLISHING.  1977.  PAGES  202.  143. 

1570  I 

1580  I  SILICON:  1 68  +  c.XP< -0 . 00458*T  ) .  T  IS  SILICON  TEMP.  IN  DEG.  C. 

1590  I  ALUMINA:  16.7  SILICON  RUBBER:  0.225  ACRYLIC:  0.208 
1600  I  HEAT  FLUX  q  FOUND  BY  STANDARD  EQUATION 
1610  I  q  -  AREA»DELTA  T/(SUM  OF  LENGTHS/K’S) 

1620  ' 

1630  I  IT  IS  ASSUMED  FOR  CALCULATIONS  THAT  THE  FOLLOWING  TEMPS.  ARE  THE  SAME 
1640  I  CHIPl  -  CHIP4  -  CHIP7 
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1650  I  CHIP3  -  CHIPS  «  CHIPS 
1660  I 

1670  I  E  IS  AREA  (SQUARE  n )  ,  F  &  S  ARE  L/K  VALUES 

1680  E-7.90E-5 

1690  F-2.42E-6 

1700  6-. 00458 

1710  S-5.73E-2 

1720  ! 

1730  Gloss ( 1  )-E*(T(61  >-T(5l  ) ) /( F*EXP( 6«T( 6 1  )  >+S  ) 

1740  Qloss(2  )-E«(T(60  )-T(52  )  )/ ( F«EXP( 6«T(  60 ) )+S ) 

1750  01os5(3)-E*<T(63  )-T( 53 >  >/< F«EXP( 6»T(  63 ) )+S ) 

1760  Qloss( 4  )-E»( T( 61  )-T ( 54  )  )/( F*EXP( 6«T( 6 1  ) )+S ) 

1770  01oss(5  >-E*(T( 62  )-T(55  )  )/ ( F«EXP( 6»T( 62 ) )+S ) 

1780  0ioss(6  )-E«( T(63  )-T(56  ) )/ ( F«EXP( G*T< 63 >  >+5  ) 

1790  Qioss(7)-E»(T(El  )-T( 57  )  )/ ( F»EXP( 6«T( 6 1  > )+S ) 

1 800  Qlossl 8  )-E«( T(64 )-T( 58  )  )/ ( F»EXP( G» T< 64 ) )+S ) 

1810  Q1o55( 9  )  =  E»( T( 63  )-T( 59  )  )/(F»EXP(G«T(  63 ) )+S ) 

1820 

1830  Qlosssurn  =  0 

1840  FOR  I»1  TO  9 

1850  Qlosssur!“01os55ur"*Clos5(  I  ; 

1860  NEXT  I 

1870  Qlo5savg=01oESSun/9 
1880  I 

1890  I  Rc  WILL  BE  CAlCU_«'TE0  ASSUMING  SILICON  .  IS  60  DEG.  C,  WHICH  IS 

1900  '  ABOUT  MID-RANGE  FOR  THE  EXPERIMENTS 

1910  I  Rc  •  <1/AREA)»SUM  OF  LENCTHS/K’S 

1920  '  ALSO  ASSUME  Ore  IS  5*4  OF  Re 

1930  RC-72S 

1940  Drc-36.3 

1950  I  Dtlid,  Dtside.  Otch;p,  Ots.  AND  Dtsink  BASED  ON  THERMOCOUPLE  CALIBRATION 

1960  Dtchip“.3 

'970  Dts-.3 

1980  Dtlid-.3 

1990  Dtside“.3 

2000  Dtsink“.3 

2010  I  SINCE  CONDUCTION  TEMP.  DIFFERENCE  •  Tchip  -  Ts.  Dt  IS  THE  RMS  OF  Dtchip 
AND  Dts 

2020  Dt-(Dtchip'2+0ts'2  )' .5 

2030  Dqloss-(01os5avQ)«((Dt/(Tflvg-Tsavo)  )“2+( Drc/Rc  )‘2  )' .5 

2040  I 

2050  I  AVERAGE  SINK  TEMPERATURE  CALCULATION 
2060  Tsink-(T( 1 1  )+T( 12  )+T( 13  )  )/3 
2070  ! 

2080  PRINT  ySIN6  ' 10X/ ’AVERAGE  SINK  TEMPERATURE  (C):  ”  .DD.DD’ iTsink 

2090  PRINT 
2100  I 

2110  1  THE  CHARACTERISTIC  LENGTH  TO  BE  USED  TO  CALCULATE  THE  NUSSELT  NUMBER 

2120  !  IS  LI  ,  WHICH  WAS  PREVIOUSLY  DEFINED 

2130  I 

2140  PRINT  USING  ’SX/’CHIP  Onet(W)  Delta  T  Nu  XUNC  IN  NU”.10A- 

2150  PRINT 
2160  I 

2170  !  CALCULATION  OF  NET  POWER.  Nu ,  Ra  AND  UNCERTAINTIES 

2180  FOR  J-1  TO  9 

2190  1  Dpower  IS  BASED  ON  ACCURACY  OF  THE  DVM  AND  THE  PRECISION  RESISTORS. 

2200  I  Ov  FOR  voltage  DROPS  ACROSS  THE  CHIP  RESISTOR  OR  THE  PRECISION  RESISTOR 
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2210  !  IS  5E-B  V.  Drp  -  0.0B  OHMS. 

2220  Dv-B.E-E 
2230  Drp-.0B 
2240  FOR  I-l  TO  3 

22  S0  Dp  oyer  {  I  )*Pouer{  I  )«((3»Dv/Enf(C>)''2+(Dv/Emf(B))“2  +  (Drp/Rp)'2)".B 

2260  NEXT  I 
2270  L-4 

2280  FOR  I-l  TO  3 

2290  Dp oyer ( L  )-Poyer( L>«((Dv/Enf<L))'2+( Dv/Enf ( I+B  )  >“2  +  ( Drp/Rp  )“2  > “ . B 

2300  L-L+1 

2310  NEXT  I 

2320  FOR  1-7  TO  9 

2330  Dp  oyer  (I  )-Poyer'  I  )»((3»0v/En?(4))''2+(  Dv/Enf  (  9  ) )  “  2  +  (  Drp/Rp  )  “2  1 '  .  B 
2340  NEXT  I 
23B0  ! 

2360  I  CALCULATION  OF  Onet 
2370  Onet( J  )-Poyer( J  )-01o9s( J  ) 

2380  Dqnet-(Dpoyer(  J  )''2+Dq1o9£'2  I"  .B 
2390  ! 

2400  !  CALCULATION  OF  Tfiln 

2410  Tf  i  ln“(  TavQ-t-Tsink  )/2 

2420  I  Dtavg  BASED  ON  TSE  CALIBRATION 

2430  Dtavg-. 27B 

2440  I 

24B0  I  CALCULATION  OF  Delta  T  BASED  ON  LID  AREA  COMPRISINS  ABOUT  BBS  OF  TOTAl. 

2460  I  CONUECTION  AREA,  AND  THE  CHIP  SIDES  4BS  OF  THE  AREA 

2470  Delt( J  )-( .BB«Tlid(J  )+ . 4B»Ts i de( J )  )-Ta ink 

2480  Ddelt-(Dtlid''2+Dtside'2+Dtsink'2  )".B 

2490  i 

2600  I  CALCULATION  OF  CONUECTION  COEFFICIENT 
2610  H( J  )-0net( J  )/( Atot»Delt( J  ) ) 

2=20  Dh<  J  )-H( J  )•( (Danet /Qnetl J  ) ) '2r < Ddel t /Del t ( J  >  > ‘2  . B 

2530  NEXT  J 
2540  I 

2560  >  PHYSICAL  PROPERTIES  ARE  TAKEN  FROM  THE  198=  3M  PRODUCT  MANUAw 

2560  I  FOR  FLUORINERT  ELECTRONIC  LIQUIDS 
2570  I 

2580  IF  D1iq-0  THEN 
2590  I 

2600  t  CALCULATION  0=  FC-75  THERMAL  CONDUCTIVITY 
2610  K-( .E5-7.89474E-4«Tf  iln  )/10 

2620  I 

2630  !  CALCULATION  OF  FC-75  DENSITY 

2640  Rho-( 1 .825-.00246»Tfiln)»1000 
2650  I 

2660  I  CALCULATION  OF  FC-75  SPECIFIC  HEAT 
2670  Cp-< .2411 1 1+3.7037E-4»Tfiln)«4187 
2680  I  THE  4187  CONVERTS  FROM  CALORIES  TO  JOULES 
2690  I 

2700  I  CALCULATION  OF  FC-75  KINEMATIC  VISCOSITY 

2710  N-1 . 4074-2. 9B4E-Z»Tfiln+3. 801 8E-4«Tfiln'2-2. 7308E-B«Tfi ln'3+8. 1B79E-9*Tfil 

n‘4 

2720  I  CONVERT  FROM  CENTISTOKES  TO  n*2/8 
2730  N-N» 1.^-6 
2740  I 

2750  I  CALCULATION  OF  THE  COEFFICIENT  OF  THERMAL 
2760  !  EXPANSION  IBETAI 
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3330  Nu< J '-H( J  )»L1/K 
3340  Dnu-Nm  J  )*(Dh(  J  )/H(  J  ) ) 

3350  Pernu«( Dnu/Nu( J )  )» 1 00 
3350  ! 

3370  !  CALCULATION  OF  6RASH0F  NUMBER 

3380  6r<J  )-9.81«BetB*<L1‘'3)»Delt(J  )/N“2 
3390  Dgr-Br’ J  )»(Ddelt/Delt( J ) ) 

3400  Pergr-C  Do''/'Sr(  J  )  )•  1  00 
34 10  ! 

3420  !  CALCULATION  OF  RAYLEIGH  NUMBER 

3430  Ra( J  )-6r( J  )*Pr«l .E-B 
3440  Dra"Ra( J )•( Dgr/6r{ J ) ) 

3450  Parra"; Dra/Ra( J ))• 1 00 
3460  I 

3470  I  DATA  AND  UNCERTAINTY  OUTPUT 

3460  IF  J-B  THEN 

3490  PRINT 

3500  PRINT 

3510  PRINT 

3520  PRINT 

3530  PRINT 

3540  PRINT 

3550  PRINT  USING  ‘lOX/’RAGE  2  OF  FILE:  "  *  ,  1  0A"  lOldf  i  leS 
35B0  PRINT 

3570  PRINT  USING  '  1 0X .D , 1 X .3(  5X ,DD . DD , > .5X ,DDD . DDD“ i J  ,Qnet ( J  )  .Del t ( J >  ,Nu( J  )  ,Per 
nu 

3580  ELSE 
3590  PRINT 

3600  PRINT  USING  *  1 0X .0  , U ,3( 5X ,DD . DO . ) ,5X .DDD . ODD' » J .Onet ; J ) ,Oel t ( J ) ,Nu( J > .Per 
nu 


3G10  END  IF 
3620  PRINT 
3G30  PRINT  USING 
J  ) 

3640  PRINT  USING 
J  ) 

3650  PRINT 
3:S0  PRINT  USING 
D . DDD ' I Dnu 
3670  PRINT  USING 
D.DDD" (Dr a 


•12X.““TEMP  BASED  RAYLEIGH  NUMBER  •  E-B  IS:  " '  .DDDDDD . DDD* t Ra( 
*12X  .""GRASHOF  NUMBER  IS:  " * .DDDDDDD.DDD* i6r( 

•6X.*“UNC  IN  THE  NUSSELT  NUMBER  IS:  **.D 

*6X.“*UNC  IN  THE  TEMP.  BASED  RAYLEIGH  NUMBER  •  E-B  IS:  * "  .0 


3G80  PRINT  USING  *BX,""UNC  IN  THE  GRASHOF  NUMBER  IS:  ’".DDD 

DDD.D" iDgr 

3690  I  ALGEBRAICALLY.  Perra  •  Pe-gr 

3700  PRINT  USING  ■6X.”ZUNC  IN  THE  RAYLEIGH  OR  GRASHOF  NUMBER  IS:  ”  ,D 

D.DDD" I  Parra 
3710  NEXT  J 
3720  ! 

3730  Ra9uri“0 


3740  Nusun"0 


3750  Onetsur“0 

3760  Deltsur"”0 

3770  6rsur"0 

3780  FOP  J"1  TO  9 

3790  Ra9un"Ra( J  )+Rasun 

3800  Nusun“Nu< J >+Nusun 

3810  0net9uw«0net( J  )+0net5un 
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3620  Del  1 5upi«De  1 1  (  J  )+De  1 1  SLiM 

3830  Gf'5un"6r (  J  j+Grsun 

3840  NEXT  J 
3850  Raovo“Rasun/9 
3880  NuevB"Nu9un/9 
3870  0netavB"Qnetsupn/9 
3880  DeltavB“Deltsum/9 
3890  6ravQ“Srsum/9 
3900  ! 

3910  FOR  J-1  TO  3 
3920  Rouira<J)-0 
3930  Rounuf  J  )“0 
3940  RowBr(J)“0 
3950  NEXT  J 
3960  FOR  J-;  TO  3 

3970  Roura( J )-(Ra( J  )+Ra<  J+3  )+Ra( J+B  )  )/3 
3980  Rownu( J  )“( Nu( J  )+Nu( J+3  )+Nu( J+6  )  )/3 
3990  Rouo+( J )“( Gr( J  )+&r( J+3  )+Sr( J+6  )  )/3 
4000  NEXT  J 
4010  I 
4020  PRINT 

4030  PRINT  USING  MBX/’TOP  ROU  FiUG  Ra«E-6  IS:  *  “  .DDDDD.DDD"  iRoura(  3  ) 
4040  PRINT  USING  •19X,”TOP  ROU  i«iU6  Nu  IS:  •  •  .DDD .  DD’ ;  Rounu(  3  ) 

4050  PRINT  USING  ’19X/‘TOP  ROW  AUG  Gr  IS:  *  “  .DDDDDDD  .  DDD“  i  Roygrl  3  > 

4060  PRINT 

4070  PRINT  USING  MBX/'MID  ROU  AUG  Ra«E-G  IS:  *  ’  .DODDD.DDD*  iRowra(  2  ) 
4080  PRINT  USING  ’lOX/’MID  ROU  AUG  Nu  IS:  ’  -  .DDD  .DD- 1  Rownu  (  2  ) 

4090  PRINT  USING  'igX.-’MIO  ROU  AUG  Gr  IS:  *  .DDDDDDD.DDD' tRougrl  2  ) 

4100  PRINT 

4110  PRINT  USING  ‘igx.’-BOT  ROW  AUG  Ra*E-B  IS:  ’ *  .DODDD.DDD' iRowral 1  ) 
4120  PRINT  USING  ‘igX.’-BOT  ROW  AUG  Nu  IS:  * "  .DDD.DD” tRownu( 1  ; 

4130  PRINT  USING  ’igX.-'BOT  ROW  AUG  Gr  IS:  " “  .DDDDDDD . DDD" ;Rouor ( 1  ) 

4140  PRINT 

4150  PRINT  USING  "igX.’-ARRAY  AUG  Onet  IS:  " “  ,D . DD’ lOnet avg 

4160  PRINT  USING  •19X,’'ARRAY  AUG  Delta  T  IS:  ■ ’  ,DD. DD* (Del tevg 

4170  PRINT  USING  "igx.’-ARRAY  AUG  Ra»E-6  IS:  *  *  .DDDDD . DDD* (Raavg 

4180  PRINT  USING  ’19X.**ARRAY  AUG  Nu  IS:  *  *  ,DDD . DD" ( Nuavg 

4190  PRINT  USING  *19X,*"ARRAY  AUG  Gr  IS:  *  *  .DDDDDDD . DDD" (Grevg 

4200  PRINT  USING  *19X,’*ARRAY  AUG  Pr  IS:  *  *  ,DDDD . D’ i Pr 

4210  ASSIGN  eFile  TO  * 

4220  END 
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APPENDIX  E.  UNCERTAINTY  ANALYSIS 


The  software  programs  CALCDIEL  and  CALC2  also  calculated  a  standard 
zeroth  order  uncertainty  analysis  as  described  in  Beckwith  and  Marangoni 
(1990)  .  The  uncertainties  were  calculated  for  the  Nusselt  number  and 
Rayleigh  number  for  each  component  on  the  circuit  board.  A  complete 
description  of  the  CALCDIEL  uncertainty  analysis  can  be  found  in  Matthews' 
thesis.  The  analysis  performed  by  CALC2  for  the  NSWC  circuit  board  is 
described  below  in  expression  form,  with  numerical  results  omitted  for 
generality.  Thi  small  uncertainty  associated  with  the  thermophysical 
properties  was  neglected.  Other  assumptions  are  included,  where 
appropr iate . 


1 . 


Conduction  heat  loss  through  the  circuit  board  assembly 


0 


J.OSS 


Arc 

Rc 


^^loss  ~  ^ 


loss' 


N 


(6RA 

1 

\  / 

AT,.  =  T... 


cat,  =  +  6r 


5AT....  =  6t„  =  ±0.3  °c 

R,  was  calculated  for  T...,,p  =  60  °C.  SR^  was  then  assumed  to  be  5%  of  R^ . 

5R,.  =  ±36.3  “C/W 

2.  Power  supplied  to  the  chip  resistors 


i  Power  =  Power 


\ 


+  — 


IP 


IP 


\  ^p  / 
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±0.000005  V 


=  5v^.  = 

8Rf,  =  ±0.05  O 

3 .  Net  power  dissipated 


Qn*t  =  Power  -  Q,,53 

=  sjfiPowei^  +  ^QIoss 


4.  Temperatures  for  calculation  of  average  heat  transfer  coefficient 

=  ±0.275  °C 
ST,.„,  =  ±0.3  °C 
AT  =  (0.55T;,^^  +  0.45T3,.^)  - 

55%  of  the  exposed  chip  surface  area  is  the  lid,  and  the  chip  sides  mak 
up  the  remainder.  Tsid<  is  based  on  the  average  between  the  li 
temperature  and  the  indicated  chip  temperature. 

8Ar  =  ± 


5  . 


dT-,.  =  5t. 

,,,  =  5t,.„,  -  ±0.3  ° 

transfer  coefficient 

h  - 

AcoAT 

6h  = 

^  I  8Ar\2 

N 

[On.c)  lAr/ 

The  small  uncertainty  in  A...  is  neglected. 
6.  Nusselt  number 


Nu  = 


hL 

k 
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a  'V 


The  small  uncertainty  in  L  is  neglected,  and  8k  is  assumed 
negligible.  Therefore 

6Nu  =  Nu^ 
h 


7.  Rayleigh  number 


Ra  =  GrPr 


V2 


Pr  = 


a 


Neglecting  property  and  length  uncertainties, 

fiGr  = 

At 

bRa  = 

Gr 


to  be 
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APPENDIX  F.  TSE  CAlilBRATION 


The  entire  NSWC  circuit  board  assembly  was  placed  inside  a  400  W 
Central  Scientific  Company  oven  for  calibration.  The  oven  had  a 
cylindrically  shaped  internal  volume  of  0.11  m’.  A  platinum  resistance 
thermometer  was  inserted  through  a  small  hole  in  the  top  of  the  oven.  A 
Rosemont  Engineering  Company  galvanometer  and  commutating  bridge,  accurate 
to  ±0.0001  ohm,  was  used  to  measure  the  thermometer's  resistance. 
Temperature  was  then  read  from  pre-printed  NBS  calibration  data,  accurate 
to  rO.Ol  °C.  Due  to  the  expected  nature  of  the  dielectric  liquid 
temperature  fluctuations  during  the  experimental  runs,  the  calibration 
data  was  rounded  to  the  nearest  0.1  °C. 

This  temperature  was  considered  to  be  the  reference  tem-.erature  for 
the  calibration  curve.  The  voltages  were  read  with  a  HP-3456A  digital 
voltmeter,  accurate  to  ±0.000005  V.  Data  points  were  taken  at  four 
different  oven  settings,  ranging  from  19.5  to  95.8  °C.  Additional  data 
points  were  taken  as  necessary  to  ensure  steady  state  conditions  had  been 
reached.  The  data  points  used  for  the  TSE  calibration  are  included  in  the 
table  below; 


channel 

60 

61 

62 

63 

64 

Temp . 

( °C; 

TSE  « 

n 

4 

5 

6 

8 

V,,(V) 

— 

0.96907 

0.96827 

0.96925 

0.96954 

0.56955 

19 . 5 

V,,(V) 

0.92845 

0.92752 

0.92861 

0.92875 

0 . 92892 

43  .  6 

V,,(V) 

0.88470 

0.88367 

0.88486 

0.88479 

0.88521 

68 . 6 

Vb.(V) 

0 . 83666 

0.83551 

0.83683 

0.83697 

0.83724 

mm 
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